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NOMENCLATURE 

A^  «  Cross-sectional  area  which  contributes  to  dynamic  inertia. 

Ag  »  Cross-sectional  area  which  contributes  resistance  to  shearing. 

Cp  *  Plate  velocity  •  {E/[v(l-v2)  J}5* 

4 

Ci  *  Dilatational  wave  velocity  in  a  beam  *  [EI^/pI^J 

Cz  *  Shear  wave  velocity  in  a  beam  =  [A  G/pA  K 

s  X 

h 

Ci  *  Shear  wave  velocity  in  a  plate  *  { G/p ] 

D  »  Flexural  rigidity  of  plate  *  Eh3/[12(l-v2) ] 

E  *  Modulus  of  elasticity 
G  ■  Modulus  of  rigidity  *  E/[2(l+v)] 
h  »  Plate  thickness 

Ib  *  Moment  of  inertia  which  contributes  resistance  to  bending 
1^  *  Moment  of  inertia  which  contributes  resistance  to  dynamic  inertia 
j  *  Superscript  referring  to  quantities  of  the  jth  cell 
k2  *  Shear  correction  factor  for  plate 
M  ■  Internal  bending  moment 
Mf  »  Radial  bending  moment  per  unit  length 
MQ  ■  Tangential  bending  moment  per  unit  length 
Qr  *  Transverse  shear  force  per  unit  length 
q  *  Intensity  of  distributed  external  load  on  beam 
r  -  Radial  distance  along  plate 


*•  - 


X 


t  *  T  ime 

V  *  Vertical  shear  force  on  a  cross-section  of  the  beam 
v  =  Velocity  of  deflection  in  a  beam  or  plate,  y  or  wt>  respectively 
w  =  Transverse  displacement  of  the  midplane  of  plate 
x  =  Coordinate  along  length  of  beam 
y  =  Deflection  of  beam 

=  Total  slope  of  the  deflection  curve  of  the  beam 

e ,  *  Angular  strain  of  an  element  of  the  beam 

*  Total  slope  of  the  deflection  curve  of  the  plate 

e,  *  Angular  strain  of  an  element  of  the  plate 
<p 

9  *  Tangential  direction 

v  *  Poisson's  ratio 

p  *  Density  of  the  material  of  beam  or  plate 

p  *  Rotation  of  the  cross-section  of  the  plate  about  the  tangential  axis 

\|>  *  Slope  of  the  deflection  curve  of  a  beam  when  shearing  force  is  neglected 

w  *  Angular  velocity  of  rotation  of  an  element  of  the  beam  or  plate, 
or  $  ,  respectively 
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1.  INTRODUCTION 


1.1  Background  Information 

The  current  Air  Force  Nondestructive  Pavement  Testing  (NDPT)  method 
contains  two  main  components — data  collection  equipment  and  analytical 
method  [1].  The  data  collection  equipment  is  used  primarily  for  evalua¬ 
tion  of  in-situ  elastic  modulus  of  each  pavement  constituent  material. 

The  elastic  modulus  values  are  then  used  as  input  data  for  the  analytical 
method  which  determines  the  structural  capacity  of  existing  airfield 
pavements . 

The  data  collection  equipment  consists  of  an  impulse  loader  with 
the  necessary  instrumentation  and  a  desk-top  computer  for  preliminary 
data  analysis  and  evaluation.  The  pavement  response  to  the  impulse 
loading  is  measured  by  using  accelerometers;  the  acceleration-time  data 
are  analyzed  by  using  the  Discrete  Fourier  Transform  (DFT)  technique  [2] 
to  obtain  the  phase  angle  versus  frequency  relationship  which  is  then 
used  to  develop  the  dispersion  curve.  From  the  dispersion  curve,  the 
shear  wave  velocity  propagating  through  the  pavement  is  obtained  and 
the  shear  modulus  of  each  pavement  constituent  material  is  computed. 
Although  the  computation  of  elastic  modulus  is  simple,  the  interpreta¬ 
tion  of  dispersion  curves  for  modulus  calculation  is  not  easy  and 
straightforward  [3,4].  For  this  reason,  the  current  NDPT  method  is 
without  an  adequate  method  for  computation  of  layer  modulus  values 
[5]. 
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Pavement  response  to  dynamic  loading  is  a  complex  problem  due  to 
the  presence  of  layers  of  different  materials  and  also  the  nature  of 
the  constituent  materials.  The  presence  of  layers  causes  reflection 
and  diffraction  of  waves.  Under  dynamic  loading,  the  pavement  materials 
such  as  bituminous  concrete  and  soil  may  produce  real  and  imaginary 
modulus  components  which  vary  with  frequency  of  vibration  [6] .  Because 
of  these  factors,  it  is  very  difficult  to  find  a  solution  especially 
for  high  order  modes  of  response  [7].  As  a  consequence,  no  straightforward 
procedure  is  presently  available  for  determination  of  realistic  modulus 
values  from  dynamic  loading  test  results  [1]. 

However,  there  is  the  so-called  method  of  Direct  Analysis  which 
bypasses  the  explicit  use  of  differential  equations  and  may  avoid  some 
of  the  difficulties  just  mentioned.  Because  this  method  has  been 
successfully  applied  to  certain  beam  and  plate  problems  [8,9],  this 
research  is  undertaken  to  investigate  the  feasibility  of  using  this 
technique  to  analyze  the  response  of  elastic  layered  system  to  dynamic 
loading.  While  the  ultimate  goal  of  the  research  is  to  develop  a 
computer  program  for  elastic  layered  system,  the  immediate  objective  of 
this  study  is  specified  below. 

1.2  Research  Objective 

The  ultimate  goal  of  the  research  was  to  develop  a  computer 
program  for  analysis  of  the  response  of  elastic  layered  system  to 
dynamic  loading.  As  a  first  step  striving  toward  the  goal,  this 
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research  was  undertaken  to  test  the  available  theory  and  to  develop  the 
modifications  necessary  for  direct  application  of  the  Method  of  Direct 
Analysis  to  problems  of  infinite  beams  and  plates  on  elastic  foundation. 


, 
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2.  INFINITE  BEAM  ON  ELASTIC  FOUNDATION 


2.1  The  Physical  Laws 

The  physical  laws  used  in  the  development  of  computer  programs 
include  the  impulse-momentum  laws  and  the  constitutive  relations.  The 
beam  is  divided  lengthwise  into  elements  of  equal  length.  For  each 
element,  the  equations  of  motion  and  constitutive  relations  are  derived 
for  two  modes  of  motion  —  rotation  and  translation.  The  basic  assumptions 
required  in  the  derivation  of  the  physical  laws  include  that  the  material 
is  elastic,  homogeneous  and  isotropic,  that  the  deformation  of  the  beam 
is  due  to  both  flexural  motion  and  shear  deflection  of  the  cross-section, 
and  that  the  elastic  foundation  behaves  as  Winkler  foundation.  These 
basic  assumptions  may  be  found  elsewhere  110]. 

2.1.1  Equations  of  Motion 

The  equations  of  motion  are  obtained  by  applying  the  impulse-momentum 
laws  directly  to  an  element  of  the  beam.  Figure  1  shows  the  free-body 
diagram  of  a  typical  jth  element.  For  this  element,  the  rates  of  change 
in  both  rotation  and  translation  are  given  below: 

Rotation  of  the  jth  element  - 

5T  '  11/2  (vi  +  v3+1)  +  ^ 


(1) 
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Translation  of  the  jth  element  - 


dv 

dt 


__L_  (vj+1  -  VJ  +  qdx  -  R.) 
pA.dx  j 


(2) 


The  foundation  reaction,  R  ^ ,  shown  in  Eq.  (2)  is  a  function  of 
foundation  deformation  which,  in  essence,  is  the  deflection  of  the  beam. 
Two  different  foundation  supports  are  considered  in  the  computer  program 
i.e.,  linear  and  nonlinear  supports.  For  the  linear  support,  the  founda 
tion  reaction  is  directly  proportional  to  the  deflection,  viz 


Rj  -  -kyj  (3) 

For  the  nonlinear  support,  the  following  hyperbolic  function  is  used: 

Rj  -  -ky^/d  +  nyj)  (4) 

where  k  is  the  spring  constant  having  a  unit  of  psi,  y^  is  the  mean 
deflection  of  the  jth  element  and  n  is  a  constant  with  a  unit  of  in 
The  above  hyperbolic  relationship  between  load  and  deformation  is  often 
used  for  describing  the  nonlinear  behavior  of  soils  [11,12]. 

2.1.2  Constitutive  Equations 

The  constitutive  equations  are  obtained  from  the  consideration  of 
the  deformation  of  a  typical  beam  element.  As  shown  in  Figure  2, 
the  jth  element  undergoes  flexural  as  well  as  shear  deformation.  Thus 
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two  equations  can  be  obtained,  one  for  shear  force  and  the  other  for 
the  bending  moment.  According  to  Timoshenko  beam  theory  [10],  the  shear 
force  is  proportional  to  the  shear  angle,  namely 

V  ■  As  G  -ij>}  (5) 

where 

ey  -  (yj+1  -  yJ)/dx  (6) 

For  the  internal  bending  moment,  the  simpler  Euler-Bernoulli  theory  provides 
the  following  relationship: 

M  .  -EIb  <s>  (7) 

where 

"  (^J+1  “  ^>/dx  (8) 

It  has  been  shown  [13[  that  any  input  to  the  beam  which  causes  it  to 
undergo  a  flexural  mode  of  motion  immediately  produces  two  different 
waves  —  dilation  and  shear  waves.  The  dilation  wave  which  is  caused 
by  discontinuities  in  w  and  M,  as  well  as  discontinuities  in  their 
higher  derivation  propagates  at  a  velocity  Cj  which  is  equal  to 


The  shear  wave  is  related  with  discontinuities  in  v  and  V,  as  well  as 
discontinuities  in  their  higher  derivatives.  The  shear  wave  velocity 
may  be  shown  to  be  [13] 


C2  = 


(10) 


2.1.3  Boundary  Conditions 

The  problem  under  investigation  is  a  beam  on  elastic  foundation  which 
has  infinite  length  and  is  subjected  to  a  vertical  dynamic  load.  The 
dynamic  load  is  distributed  over  a  small  area  which  can  be  treated  as 
a  concentrated  load.  In  the  analysis,  the  beam  is  split  into  two  at  the 
loading  point  so  that  each  semi-infinity  beam  carries  a  shear  force 
which  is  equal  to  one-half  of  the  applied  load.  Meanwhile,  for  the 
continuity  requirement,  the  slope  of  the  deformed  beam  at  the  loading 
point  is  maintained  at  a  value  of  zero. 

Furthermore,  the  shape  of  the  deflected  beam  dictates  that  deflections 
beyond  distance  L  (measured  from  the  load)  become  very  small.  According 
to  Timoshenko  [14],  the  distance  L  for  static  loading  is  approximately 
equal  to  5.5/B,  where  S  ■  4  -y k/4EI,  k  *  spring  constant  of  the  elastic 
foundation,  E  «  Young's  modulus,  and  I  =  moment  of  inertia  of  the  cross 
section  of  the  beam  with  respect  to  z  axis.  On  this  basis,  the  semi¬ 
infinite  beam  can  be  approximated  by  a  beam  with  a  finite  length  L.  Thus, 
for  the  problem  under  consideration,  the  shear  input  and  boundary  con¬ 


ditions  are 
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(11) 

(12) 

(13) 


In  the  analysis,  it  is  convenient  to  express  the  boundary  conditions  in 
terms  of  the  angular  and  linear  velocities,  respectively,  as  follows. 


(14) 

(15) 


where  subscript  t  denotes  time. 

The  preceding  equations  are  a  complete  statement  of  the  problems  of 
the  flexural  traveling  stress  waves  in  an  infinite  beam  on  elastic  founda¬ 
tions.  Of  these  equations,  the  rotational  impulse-momentum  law,  Eq.  (1), 
indicates  that  the  angular  velocity  of  rotation,  w,  is  a  function  of  both 
the  moment  and  the  shear.  Since  the  moment  is  associated  with  the  dilation 
wave  and  the  shear  with  the  shear  wave,  Eq.  (1)  exhibits  an  interaction 
of  two  waves  and  by  definition  the  angular  velocity  is  a  coupled  quantity. 
Moreover,  Eq.  (5)  indicates  that  the  shear  force  is  also  a  coupled 
quantity  since  it  is  a  function  of  both  the  linear  displacement  gradient 
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and  the  angle  of  rotation.  Because  of  the  interaction  or  coupling  between 
the  two  waves,  variables  which  are  associated  with  both  wave  speeds  must 
be  given  special  treatment.  Thus,  the  solution  requires  a  special 
technique  of  coupling  which  is  presented  in  the  following. 

2.2  The  Technique  of  Coupling 

In  the  Method  of  Direct  Analysis,  the  beam  is  divided  into  finite 
number  of  elements  or  cells  as  mentioned  earlier.  It  is  convenient  to 
determine  the  cell  length  based  on  the  "characteristic  assumption," 
namely,  dx  -  c  dt  [8,15],  where  c  *  wave  velocity  and  dt  *  time  increment. 
Since  two  distinct  wave  speeds  are  present,  the  cell  length  dx  -  c  dt 
cannot  be  satisfied  for  both  speeds  simultaneously.  Thus,  a  technique 
is  needed  to  determine  whether  the  wave  under  consideration  has  passed 
across  the  cell.  If  the  wave  has  crossed  the  cell,  the  computation 
of  stresses,  strains  and  velocities  for  this  cell  will  be  triggered. 
Otherwise,  further  computation  is  skipped  and  the  current  dynamic 
variables  are  held  over  until  the  time  has  sufficiently  advanced,  and 
the  cycle  is  repeated. 

There  are  two  procedures  available  for  carrying  out  the  above 
technique,  namely,  the  "wave  index"  and  "clock"  methods  [16].  These  two 
procedures  are  completely  analogous  but  the  "clock"  method  is  more 
amenable  to  a  problem  in  which  two  waves  are  present.  For  this  reason, 
the  "clock"  method  is  adopted  in  this  analysis  and  is  discussed  below. 


v--  TJ>V~ 
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In  the  "clock"  method,  two  time  "clocks"  are  established;  one  for 
dilatation  wave  and  the  other  for  shear  wave.  Each  "clock"  signals  the 
instant  when  the  wave  propagation  procedure  for  that  particular  wave 
is  undergone.  The  rate  of  advancement  of  the  two  clocks  are  defined 
as  follows: 


dtD 


dx 


C 


D 


for  dilatation  wave  clock,  and 


(16) 


for  shear  wave  clock, 

in  which  Cp  and  Cg  are  the  dilatation  and  shear  wave  velocities,  respectively. 

Since  Cp  is  greater  than  Cg,  the  dilatation  wave  front  initially  will 
reach  the  end  of  a  cell  before  the  corresponding  shear  wave  front.  When 
this  has  happened,  the  dilatation  time  clock  will  advance  dtp,  viz, 
tp'  *  tp  +  dtp,  and  the  quantities  associated  with  the  dilatation  wave 
will  be  adjusted.  Then,  two  possibilities  follow.  One  possibility 
is  that  the  dilatation  clock  has  advanced  another  dtp  before  the  shear 
clock.  When  this  condition  occurs,  the  quantities  associated  with  the 
dilation  wave  front  again  will  be  adjusted.  The  other  possibility  is 
that  the  shear  wave  clock  advances  dtR  before  the  next  advancement  of 
the  dilatation  clock.  When  this  happens,  the  quantities  associated  with 
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the  shear  wave  front  will  be  adjusted.  A  flow  diagram  summarizing  this 
propagation  procedure  is  given  in  Figure  3.  Note  that  in  this  analysis, 
the  actual  time,  t,  is  chosen  to  coincide  with  t^.  A  more  detailed 
discussion  on  wave  propagation  procedures  is  available  elsewhere  [16]. 

The  preceding  wave  propagation  technique  is  used  as  the  basis  for 
the  arrangement  of  the  various  physical  laws  presented  earlier.  It 
should  be  emphasized  that  the  position  of  the  important  steps  must  be 
carefully  organized  in  order  to  maintain  a  stable  solution.  In  general, 
the  arrangement  requires  that  the  cause  and  effect  are  properly  sequenced 
and  the  statement  of  boundary  conditions  should  follow  immediately  after 
the  calculation  of  the  dynamic  quantity  for  which  they  are  prescribed. 


2.3  The  Technique  of  Considering  Damping 

The  effect  of  damping  is  taken  into  consideration  with  the  use  of 
exponential  damping  function.  In  the  problem  of  flexural  traveling  waves 
in  beams  and  plates  there  are  two  velocities  —  linear  and  angular  velocities 
which  must  be  damped  in  order  to  attain  an  equilibrium  state.  According 
to  exponential  damping  function  which  is  often  used,  these  two  velocities 
can  be  expressed  as  follows: 


u> 


oj  e 
o 


-t/rl 


(18a) 


v 


v  e 
o 


-t/x 


2 


(18b) 
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Figure  3.  Flow  Diagram  for  "Clock"  Method  of  Wave  Propagation 


where  and  vq  are  the  initial  linear  and  angular  velocities,  respectively, 
and  t :  and  x,  are  arbitrary  constants. 

Upon  differentiation,  one  obtains  the  magnitude  of  damping  occurring 
in  time  dt  or  the  change  of  angular  and  linear  velocities  due  to  damping 
as  follows: 

du  =  -  w  —  (19a) 

x , 


These  damping  relations  Eq.  (19a  and  b)  are  incorporated  into  the  previously 
developed  impulse-momentum  laws  as  shown  in  Figure  4.  It  should  be  noted 
that  because  of  the  desire  to  eliminate  computation  fluctuations,  it 
has  been  found  necessary  to  use  half  of  the  initial  and  half  of  the  final 
velocity  (linear  or  angular)  to  obtain  the  total  contribution  of  motion 
damping  occurring  in  time  dt. 

From  Eq.  (19)  it  is  seen  that  a  completely  undamped  solution  can 
be  obtained  by  letting  x ^  and  x  equal  to  infinity.  To  obtain  a  transient 
and  static  solutions,  dt/x^  must  be  restricted  to  lie  in  the  range  between 
zero  and  one.  This  general  principle  of  synthesizing  static  and  dynamic 
solutions  constitutes  a  unique  feature  of  the  Method  of  Direct  Analysis. 

Note  that  in  mathematical  approaches,  static  and  dynamic  solutions  are 
governed  by  elliptic  and  hyperbolic  equations,  respectively.  Separate 
analyses  are  generally  required  for  each  type  of  problem. 
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Using  these  properly  organized  physical  laws,  boundary  conditions 
and  damping  equations,  the  computer  program  developed  by  Koenig  [16] 
for  cantilever  beams  has  been  modified  and  extended  to  suit  the  condition 
of  infinite  beams  on  elastic  foundations.  The  computer  program  which  is 
developed  for  I3M  370/3081  is  capable  of  providing  the  response  of  beams 
to  both  transient  and  steady-state  loadings.  The  transient  loading  may 
be  a  triangular,  rectangular  or  haver  sine  pulse;  and  the  steady-state 
loading  may  include  harmonic  or  sustained  loadings.  Furthermore,  both 
linear  and  non-linear  elastic  foundation  can  be  considered. 

The  computed  response  can  be  presented  in  both  numerical  and  graphical 
forms.  The  plotting  of  output  is  accomplished  by  incorporating  an 
available  program  PLOTIT  [17]  into  the  computer  program.  The  listing 
of  the  complete  computer  program  is  included  in  Appendix  A. 


3.  DYNAMIC  RESPONSE  OF  BEAM  ON  ELASTIC  FOUNDATION 


The  developed  computer  program  was  used  to  analyze  the  response  of 
infinite  beam  on  elastic  foundation  for  several  types  of  dynamic  loading 
to  demonstrate  its  effectiveness  in  the  analysis  of  stress  wave  propa¬ 
gation.  The  dynamic  loadings  analyzed  included  step  sustained  loading, 
step  pulse  loading  and  sinusoidal  loading.  Results  of  the  analyses 
are  presented  and  discussed  below. 

3.1  Response  to  Step  Sustained  Loading 

In  this  analysis,  the  infinite  beam  is  made  of  portland  cement 
concrete  having  a  cross  section  1  in.  wide  and  4  in.  high  with  a 
modulus  of  elasticity  of  3  x  10^  psi.  The  beam  is  subjected  to  a  con¬ 
centrated  vertical  load  of  200  pounds  with  a  sufficiently  long  duration 
so  that  it  can  be  considered  as  a  sustained  load.  The  beam  is  supported 
by  springs  having  two  different  stiffnesses,  one  lower  and  the  other 
higher  than  the  stiffness  of  the  beam.  To  demonstrate  the  ability 
of  the  computer  program  for  handling  different  damping  conditions, 
both  under  and  over-damped  conditions  are  analyzed  for  one  spring  stiffness. 

3.1.1  Over-Damped  System 

In  the  preceding  chapter,  it  has  been  pointed  out  that  for  over¬ 
damped  condition,  the  constant  in  Equation  (19)  must  be  small. 

It  should  be  noted,  however,  that  too  small  values  of  may  over- 
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depress  the  response  and  as  a  result  considerably  longer  time  is  needed 
to  reach  a  steady-state  condition.  The  time  required  to  reach  a  steady- 
state  response  depends  greatly  on  the  stiffness  of  the  beam-foundation 
system  and  time  increment  (dt)  used  in  the  analysis.  Thus,  the  values 
of  must  be  carefully  chosen  in  order  to  keep  the  time  to  steady-state, 
and  therefore  the  computation  cost,  to  a  minimum. 

4  6 

In  the  analysis,  two  spring  stiffnesses,  1  x  10  and  4  x  10  psi, 

4 

are  used.  For  the  system  with  1  x  10  psi  spring  constant,  the  length 
of  the  beam  analyzed  is  48  in.,  cell  length  is  1  in.  and  the  time 
constant  is  chosen  at  30  usee.  Results  of  the  analysis  are  presented 
in  two  general  forms:  response  versus  time  and  response  versus  distance. 
Figures  5,  6  and  7  show  the  deflections,  moments  and  shears  at  three 
locations  versus  time,  respectively.  It  is  seen  that  the  responses 
almost  reach  a  maximum  at  the  end  of  the  curves.  Also,  time  lags  exist 
between  the  responses  at  the  loading  point  and  that  away  from  it  indicating 
the  phenomenon  of  stress  wave  propagation  away  from  the  loading  point. 

The  computed  responses  shown  at  the  end  of  the  curves  are  plotted 
against  distance  in  Figures  8  (deflection),  9  (slope),  10  (moment) 
and  11  (shear) .  Also  included  in  these  figures  are  an  available  exact 
solution  for  static  loading  which  is  given  by  Timoshenko  and  is  described 
below  [14]. 

Deflection  y  -  ^  e  (cos  Bx  +  sin  Bx)  (20a) 


Variation  of  Beam  Deflections  with  Time  for  Three^Locations 
Measured  from  Loading  Point,  Spring  Constant  =  10  psi. 
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Figure  8.  Variation  of  Deflection  with  Distance  and  Comparison  with  Exact 
Solution. 


STEP  LOADING,  P  =  200  LB 
LOAD  DURATION,  T  =  3  SEC 
SPRING  CONSTANT,  K  =  10,000  PSI 
TIME  OF  COMPARISON,  TC  =  7.83  MSEC 


Figure  10.  Variation  of  Moment  with  Distance  and  Comparison  with  Exact  Solution 


DISTANCE  X  IN  INCHES 
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Slope  ^  e  sin  Sx  (20b) 

Moment  M  =  -  ^  e  ^X  (sin  Sx  -  cos  Sx)  (20c) 

Shear  V  *  -  £  e  cos  Sx  (20d) 

where  notations  P,  k  and  S  have  been  defined  in  the  preceding  chapter. 

It  is  obvious  that  the  exact  solution  provides  continuous  curves. 

For  clarity,  however,  only  some  values  are  selected  arbitrarily  for 
comparison.  Although  there  is  some  discrepancy  between  the  two  sets  of 
data,  the  agreement  generally  speaking  is  very  good.  Primary  causes 
for  the  difference  may  include  the  following:  (1)  At  the  time  of 
comparison,  the  steady-state  response  has  not  yet  established  as  revealed 
by  Figures  5,  6  and  7  in  which  the  response  has  not  quite  reached  the 
constant  value,  (2)  The  cell  length  of  1  in.  is  too  large.  Of  the  four 
figures  under  consideration.  Figure  9  provides  a  better  view  of  the 
effect  of  coarse  cell;  in  this  figure,  1-in.  line  segments  are  clearly 
shown.  It  is  believed  that  with  the  use  of  smaller  cell  length  together 
with  larger  time  of  comparison,  the  results  of  Direct  Analysis  should 
match  with  the  exact  solution  very  well.  The  excellent  agreement  between 
the  two  sets  of  results  indicates  that  the  developed  computer  program 
is  capable  of  providing  accurate  solution  for  infinite  beam'  on  elastic 
foundations . 

Another  system  analyzed  has  a  spring  constant  of  4  x  10^  psi. 

For  this  condition,  the  length  of  the  beam  used  is  10  in.,  cell  length 
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is  0.05  in.,  and  time  constant  is  chosen  at  6.6  usee.  Results  of 
the  analysis  are  summarized  in  Figures  12,  13  and  14  which  show  the 
variation  of  deflection,  moment  and  shear  with  time,  respectively.  As 
expected,  the  shape  of  the  curves  resembles  that  of  the  previous  system. 
Primary  reasons  for  analyzing  a  system  with  such  a  high  spring  constant 
are  two-fold:  (1)  to  demonstrate  that  the  computer  program  can  be  used 
to  analyze  the  dynamic  response  of  an  electric  beam  supported  by  a 
foundation  which  is  stiffer  than  the  beam,  and  (2)  to  provide  response 
data  of  over-damped  condition  for  comparison  with  the  data  of  under¬ 
damped  condition  which  is  presented  later. 

3.1.2  Under-Damped  System 

In  order  to  obtain  at  least  one  full  cycle  of  fluctuation  and  at 
the  same  time  to  keep  the  computation  cost  to  a  minimum,  a  very  high 
value  of  spring  constant  (4  x  10^  psi)  was  selected  for  analysis. 

This  value  of  spring  constant  was  also  used  in  the  preceding  over-damped 
condition  so  that  a  comparison  of  response  data  between  over-damped  and 
under-damped  systems  can  be  made.  As  before,  the  length  of  the  beam 
used  in  the  analysis  is  10  in.,  and  the  cell  length  is  0.05  in.  but 
time  constant  is  selected  at  30  msec.  Figures  15,  16  and  17  present 
respectively,  the  variations  of  deflection,  moment  and  shear  with  time 
for  three  different  locations.  Although  only  a  little  more  than  one 
cycle  of  data  are  is  obtained,  a  trend  that  the  response  data  fluctuate 
around  the  maximum  value  of  the  over-damped  solution  (Figures  12  through 


Variation  of  Beam  Deflection  with  Time  for  Three  Locations  Measured 
from  Loading  Point,  Spring  Constant  =  4  x  10^  psi.  Over  Damped. 


Figure  15.  Variation  of  Beam  Deflection  with  Time  for  Three  Locations 
Measured  from  Loading  Point,  Spring  Constant  =  4  x  10  psi 
Under  Damped. 


Figure  16.  Variation  of  Beam  Moment  with  Time  for  Three  Locations  Measured 
from  Loading  Point,  Spring  Constant  =  4  x  10^  psi.  Under  Damped. 
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14)  is  revealed.  Also,  because  of  the  damping,  the  peak  response  value 
decreases  with  increasing  cycle  of  fluctuation. 

3.2  Response  to  Step  Pulse  Loading 

The  step  pulse  loading  analyzed  is  a  200  lb.  vertical  concentrated 
load  having  a  load  duration  of  3  msec.  The  infinite  beam  is  of  the  same 
material  and  the  same  dimension  as  that  used  in  the  previous  condition, 
namely,  portland  cement  concrete  beam,  1  in.  wide  by  4  in.  high  with 
3  x  10^  psi  modulus.  However,  the  support  is  a  nonlinear  elastic 
foundation  which  behaves  according  to  the  hyperbolic  function  described 

by  Eq.  (4).  In  Eq.  (4),  the  values  of  constants  k  and  n  are  chosen  to 

-1  4 

be  10,000  psi  and  100  in.  ,  respectively.  With  the  use  of  1  x  10  psi 

for  k  value  and  a  load  of  200  lb.,  this  system  is  different  from  the 

previous  one  only  in  load  duration  and  nonlinear  spring  support.  To 

obtain  an  over-damped  solution,  the  same  time  constant  t  =  30  usee  as 

used  before  is  adopted.  Meanwhile,  the  length  of  the  beam  analyzed  is 

also  48  in.  and  the  cell  length  is  1  in. 

The  deflection,  moment  and  shear  versus  time  relationships  obtained 
from  the  computer  analysis  are  presented  in  Figures  18,  19  and  20, 
respectively.  It  is  seen  that  during  the  loading  period  the  response 
increases  with  time  following  the  same  path  as  that  of  the  step  sustained 
loading  shown  in  Figures  5,  6  and  7.  After  the  load  is  removed,  i.e., 
after  3  msec,  the  response  gradually  diminishes  with  increasing  time. 

A  comparison  between  Figures  4  through  7  and  Figures  18  through  20 


Figure  18.  Variation  of  Beam  Deflection  vith  Time  for  Three  Locations  Measured 
from  Loading  Point,  Nonlinear  Spring  Support. 


Figure  19.  Variation  of  Moment  with  Time  for  Three  Locations  Measured  from 
Loading  Point,  Nonlinear  Spring  Support. 


Figure  20.  Variation  of  Shear  with  Time  for  Three  Locations  Measured  from 
Loading  Point,  Nonlinear  Spring  Support. 
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indicates  no  apparent  difference  in  the  response  within  the  loading 
period  between  the  linear  and  nonlinear  systems.  This  is  probably  because 
the  magnitude  of  deflection  under  200  lb.  is  not  large  enough  to  cause 
significant  difference  in  foundation  reaction. 

To  show  the  variation  of  response  with  distance,  a  time  of  3.48  msec 
is  used  for  plotting.  The  curves  thus  obtained  are  presented  in  Figures  31 
(deflection) ,  22  (slope) ,  23  (moment)  and  24  (shear) .  The  shapes  of 
these  curves  in  general  are  similar  to  those  in  Figures  8  through  11 
except  for  the  shear  curve  especially  in  the  region  near  the  loading 
point.  At  x=o,  the  shear  equals  100  lb.  in  Figure  11  and  0  lb.  in  Figure  24. 
This  is  due  to  the  fact  that  the  time  used  for  plotting  (3.48  msec)  exceeds 
the  duration  of  the  pulse  loading  (3  msec).  At  3.48  msec,  the  loading 
has  already  vanished  and  therefore  the  shear  at  x=o  equals  zero  as  indicated 
in  Figure  20. 

3.3  Response  to  Sinusoidal  Loading 

The  same  portland  cement  concrete  beam,  1  in.  wide  by  4  in.  high, 
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supported  by  a  linear  elastic  foundation  with  k  *  1  x  10  psi  is  subjected 
to  a  sinusoidal  loading  which  has  an  intensity  of  200  lb.  with  a  frequency 
of  excitation  of  3600  t.  The  duration  of  the  sinusoidal  loading  is  3  sec. 
which  is  large  enough  to  provide  a  steady-state  response.  As  before, 
the  length  of  the  beam  analyzed  is  48  in.  and  the  time  constant  is 
selected  at  30  psec  to  effect  over-damping.  However,  to  reduce  the  cost 
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Variation  of  Beam  Deflection  with  Distance,  Nonlinear  Spring 
Support . 


STEP  LOADING,  P  =  200  LB 

LOAD  DURATION,  T  =  3  MSEC 

TIME  FOR  PLOTTING,  TP  =  3.48  MSEC 
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Figure  22.  Variation  of  Slope  with  Distance,  Nonlinear  Spring  Support. 


Figure  23.  Variation  of  Moment  with  Distance,  Nonlinear  Spring  Support. 


Figure  24.  Variation  of  Shear  with  Distance,  Nonlinear  Spring  Support. 
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of  computation  a  cell  length  of  2  in.  instead  of  1  in.  is  used  in  this 
analysis.  The  computed  deflection,  moment  and  shear  are  shown  respectively 
in  Figures  25,  26  and  27  as  a  function  of  time.  As  would  be  expected, 
delayed  responses  at  points  away  from  the  load  are  shown. 

The  responses  obtained  for  three  different  times,  8.7,  17.4  and 
26.1  usee  are  plotted  against  distance  in  Figures  28  through  31.  The 
time  used  for  plotting  is  chosen  arbitrarily.  These  figures  demonstrate 
again  the  phenomenon  of  wave  propagation. 
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Figure  25.  Variation  of  Beam  Deflection  with  Time  for  Three  Locations  Measured 
from  Loading  Point. 


STEP  LOADING 
P  =  200  Sin  (3600t) 
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Figure  26.  Variation  of  Beam  Moment  with  Time  for  Three  Locations  Measured 
from  Loading  Point. 


STEP  LOADING 
P  =  200  Sin  (3600t) 
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Figure  27.  Variation  of  Beam  Shear  with  Time  for  Three  Locations  Measured  from 
Loading  Point. 


STEP  LOADING,  P  =  200  Sin  (3600t) 
LOAD  DURATION.  T  =  3  SEC 
SPRING  CONSTANT,  K  =  10000  PS1 


47 


Figure  28.  Variation  of  Deflection  with  Distance  for  Three  Points  of  Time  After 
Load  Application. 
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Figure  30.  Variation  of  Moment  with  Distance  for  Three  Points  of  Time  After 
Load  Application. 


Figure  31.  Variation  of  Shear  with  Distance  for  Three  Points  of  Time  After 
Load  Application. 
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4.  INFINITE  PLATE  ON  ELASTIC  FOUNDATION 


4.1  Development  of  Computer  Program 

The  problem  of  infinite  plate  on  elastic  foundation  with  vertical 
loading  is  axisymmetrical.  For  this  reason,  in  the  development  of  the 
computer  program,  the  plate  is  divided  into  cells  or  elements  having 
the  shape  of  concentric  rings.  As  for  the  beam,  the  equations  of  motion 
and  constitutive  equations  of  each  cell  are  derived  for  both  rotation 
and  translation  motions.  Also,  as  before,  the  elastic  foundation  is 
assumed  to  behave  as  Winker  medium. 


In  the  development  of  physical  laws,  we  consider  the  free-body 
diagram  of  a  typical  jth  element  shown  in  Figure  32  and  apply  the 
impulse-momentum  laws.  The  equations  of  motion  thus  obtained  are  as 
follows: 


Rotation  of  the  jth  element 


du 

dt 


M  J+1(rj+dr)  -  M  jrJ  -  dr  -  ^  [r^  Q  j  dr  +  (rj+dr)  Q  j+1dr]  (21) 
r _ r _ 8  L _ r _ t 

(0h3/12)[rjdr  +  (dr)2/2] 


Translation  of  the  jth  element 


,  (rj  +  dr)  Q  J+1  -  rjQ  J  -  (rJ  +  dr/2)  R^ 
av  _ r _ r _ 

dt  “ 


ph  [r^  dr  +  (dr)2/2] 


(22) 


where  the  foundation  reaction  R^  has  been  defined  in  Equations  (3)  and  (4), 


Diagram  of  Plate  Element  (typical  j  element) 
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The  constitutive  equations  are  obtained  in  a  manner  similar  to 
those  of  the  beam,  viz 


M  -  D{e  +  £ 
r  or 

(23) 

Mq  *  D{—  <J>  +  v  e  } 
ur  $ 

(24) 

Qr  =  K22  Gh  4  +  ew} 

(25) 

Eh3 

0  =  12  (1-V) 

(26) 

-  (4^+1  -  <f>^)/dr 

(27) 

e  ”  (w^+1  -  w^)/dr 
w 

(28) 

K2a  =  0.76  +  0.3  v 

(29) 

As  for  the  beam,  the  dilatation  wave  velocity  is  the  velocity  of 
propagation  of  discontinuities  in  the  moments  as  well  as  the  angular 
velocity  and/or  their  higher  derivatives,  and  the  shear  wave  velocity 
is  the  velocity  of  propagation  of  discontinuities  in  the  shear  and  linear 
velocity  and/or  their  higher  derivatives.  These  two  velocities  can  be 
expressed  as  follows: 


Dilatation  wave  velocity  or  plate  velocity 


C 

P 


E 

P(l-v2) 


4 


(30) 
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Shear  wave  velocity 

K2C2X  '  K2  [  f  (31) 

The  boundary  conditions  for  the  infinite  plate  on  elastic  foundation 
may  be  derived  from  the  shape  of  the  deflected  plate.  Under  a  vertical 
load  uniformly  distributed  over  a  small  area,  the  deformed  plate  resembles 
a  bowl  which  is  axisymmetric  about  the  loading  axis  so  that  at  the  loading 
center  the  slope  equals  zero.  Furthermore,  according  to  Alpan  and 
Leshchinsky  [18],  the  plate  deflections  beyond  a  distance  of  R  measured 
from  the  loading  center  are  very  small  and  can  be  neglected  for  practical 
purposes.  The  distance  R  is  a  function  of  the  radius  of  the  loaded  area 
and  the  stiffness  of  the  system  and  can  be  estimated  from  the  following 
equation: 


R  =  Hf(r  /*) 
o 


(32) 


where 


r  *  radius  of  loaded  area 
o 

l  *  characteristic  length 


Other  notations  have  been  defined  earlier.  Also,  available  solutions  [19,20] 
indicate  that  the  deflected  shape  of  the  plate  is  oscillatory  and  has 
its  first  two  zeros  at  3.92  l  and  8.36  l.  In  this  analysis,  the  larger 
value  of  8.36  l  and  that  computed  from  Gq.  (32)  is  used  to  determine  the 
radius  of  the  plate. 


With  the  use  of  this  distance,  the  infinite  plate  may  be  approximated 
by  a  circular  plate  having  a  radius  equal  to  R  and  the  loading  acting 
at  the  center  of  the  plate.  Thus,  for  the  problem  under  consideration, 
the  shear  input  and  boundary  conditions  are 


Deflection  y  ]  T  =  v]  T  =  0 
x=“L  x=L 


(36) 


where  rQ  is  the  radius  of  the  loaded  area,  w  and  v  are  angular  and 
linear  velocities,  respectively,  and  t  denotes  time. 

Using  the  preceding  physical  laws,  constitutive  equations,  boundary 
conditions  and  the  damping  equations  described  in  an  earlier  chapter,  the 
computer  program  developed  by  Koenig  [16]  for  circular  plates  with  clamped 
edges  are  modified  and  extended  to  suit  the  condition  of  infinite  plates 
on  elastic  foundation.  The  final  computer  program  is  included  in  Appendix 
B. 


4.2  Response  to  Step  Sustained  Loading 


The  developed  computer  program  is  used  to  analyze  the  response  of 
a  4-in.  thick  infinite  plate  on  elastic  foundation  to  a  step  pulse  loading, 
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the  loading  has  a  magnitude  of  100  lb.  and  is  uniformly  distributed 
over  a  circular  area  0.2  in.  in  radius.  The  duration  of  the  loading  is 
4  sec.  which  is  long  enough  to  be  considered  as  a  sustained  loading.  The 
plate  has  a  Young's  modulus  of  3  x  10^  psi  and  a  Poisson's  ratio  of  0.20. 
The  supporting  foundation  is  a  weightless  linear  spring  having  a  spring 

4 

constant  of  1  x  10  pci.  In  the  analysis,  the  4-in.  thick  infinite  plate 
is  approximated  by  a  circular  plate  with  a  radius  of  50  in.  A  cell  length 
of  0.2  in.  is  used  and  a  time  constant  of  30  usee  is  selected  in  order 
to  obtain  an  over-damped  response. 

The  analyzed  deflection,  moment  and  shear  are  plotted  against  time 
in  Figures  33,  34  and  35,  respectively.  As  would  be  expected,  the  response 
at  a  point  increases  with  increasing  time.  In  general,  the  shape  of  the 
curves  resembles  that  of  the  curves  for  beam  on  elastic  foundation. 

4.3  Response  to  Sinusoidal  Loading 

An  analysis  is  also  made  for  the  response  of  the  same  plate  to  a 
sinusoidal  vertical  loading.  The  intensity  of  the  loading  is  also  100  lb., 
the  frequency  of  excitation  is  36000  t.  Other  conditions  including  plate 
modulus,  Poisson's  ratio,  spring  constant,  cell  length  and  time  constant 
are  the  same  as  that  used  for  the  beam  analyzed  previously.  Results 
of  the  analysis  are  presented  in  Figures  36  (plate  deflection),  37  (moment) 
and  38  (shear) . 


LOCATIONS 


vT 


Loading  Point,  Spring  Constant 


Figure  34.  Variation  of  Plate  Moment  with  Time  for  Three  Locations  Measured 
from  Loading  Point,  Spring  Constant  =  10^  psi. 


Figure  35.  Variation  of  Plate  Shear  with  Time  for  Three  Locations  Measured 
from  Loading  Point,  Spring  Constant  =  10^  psi. 


Figure  36.  Variation  of  Plate  Deflection  with  Time  for  Three  Locations 
Measured  from  Loading  Point,  Spring  Constant  =  104  psi. 


Loading  Point 


Variation  of  Plate  Shear  with  Time  for  Three  Locations  Measured 
from  Loading  Point,  Spring  Constant  =  10^  psi. 
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5.  SUMMARY  AND  CONCLUSIONS 


As  a  step  toward  the  development  of  a  computer  program  for  analyzing 
the  response  of  elastic  layered  systems  to  dynamic  loading,  the  Method 
of  Direct  Analysis  was  used  to  analyze  the  dynamic  response  of  infinite 
beams  and  plates  on  elastic  foundation.  The  Method  of  Direct  Analysis 
uses  only  the  impulse-momentum  laws  and  constitutive  relations  but  by¬ 
passes  the  explicit  use  of  differential  equations.  As  a  result  of  the 
research,  two  computer  programs  were  developed,  one  for  infinite  beams 
and  the  other  for  infinite  plates  both  on  elastic  foundations.  The 
computer  programs  are  capable  of  handling  various  types  of  dynamic 
loading  such  as  sinusoidal,  impulse,  and  sustained  loadings.  Also, 
both  linear  and  nonlinear  spring  supports  can  be  considered. 

To  demonstrate  the  effectiveness  of  the  developed  computer  programs, 
response  analyses  were  made  for  sustained,  pulse  and  sinusoidal  loadings. 
For  the  sustained  loading  on  an  infinite  beam,  both  over-damped  and 
under-damped  conditions  were  analyzed.  Also,  under  the  pulse  loading, 
the  beam  was  analyzed  for  nonlinear  spring  support.  The  results  for 
sustained  loading  on  an  infinite  beam  were  compared  with  an  available 
exact  solution.  An  excellent  agreement  with  the  exact  solution  was 
obtained. 

The  results  of  the  study  indicate  that  the  Method  of  Direct  Analysis 
is  an  effective  tool  for  dynamic  response  analysis.  It  appears  feasible 
that  a  computer  program  for  response  analysis  of  elastic  layered  systems 
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to  dynamic  loading  can  be  developed  based  on  the  principle  of  the  Method 
of  Direct  Analysis. 
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6.  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 


The  computer  programs  developed  from  this  research  are  based  on 
the  assumption  that  the  elastic  foundation  is  a  weightless  Winkler  medium. 
Therefore,  further  research  is  needed  to  extend  and  modify  the  computer 
programs  to  suit  elastic  layered  systems.  To  achieve  the  ultimate  goal 
of  developing  a  program  for  nondestructive  pavement  testing  (NDPT) ,  the 
following  steps  of  research  are  recommended: 

1.  The  dynamic  response  of  a  composite  plate  consisting  of  two 
different  materials  and  supported  by  a  Winkler  medium  should 
be  studied.  Results  of  the  analytical  investigation  should 
be  validated  experimentally. 

2.  Results  of  the  study  on  composite  plate  should  be  extended 

to  the  condition  of  an  infinite  plate  supported  by  an  Isotropic 
elastic  half-space.  Field  testing  should  be  conducted  to 
validate  the  computer  program  thus  developed. 

3.  Extend  the  computer  program  developed  above  to  the  condition 
of  a  composite  plate  on  an  isotropic  elastic  half-space.  This 
phase  of  research  should  also  be  accompanied  by  full-scale 


field  testing. 
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APPENDIX  A 

Computer  Program  for  Infinite  Beam 
on  Elastic  Foundation 


//  EXEC  PGM-IEFBR14 

//D  DD  DSN-MEN. P11850.A3B. PLOT, DISP-(OLD, DELETE) , 

//  V0L-REF-MEN.P11B50.A3B.LIB 

//D  DD  DSN-MEN. P11850. A3B .PL0T1 , DISP-(OLD, DELETE), 

//  V  OL-REF-MEN. PI 1 850 . A3B . LIB 

//D  DD  DSN-MEN. PI 1850. A3B . PL0T2 , DISP-( OLD .DELETE) , 

//  V0L-REF-MEN.P11850.A3B. LIB 

//D  DD  DSN-MEN. P11850.A3B.PLOT3,DISP-(OLD, DELETE), 

//  V OL-REF-MEN. Pll 850. A3B. LIB 

//D  DD  DSN-MEN. PI 1850. A3B.PLOT4,DISP-(OLD, DELETE) , 

//  V OL-REF-MEN. Pll 850. A3B. LIB 

//D  DD  DS N-MEN. Pll 850. A3B.PLOT5,DISP-( OLD, DELETE) , 

//  V OL-REF-MEN. P11850. A3B . LIB 

//D  DD  DSN-MEN. P11850.A3B.PL0T6,DISP-(0LD, DELETE) , 

//  V0L-REF-MEN.P11850.A3B.LIB 
//  EXEC  FGCLG, PARM-'NOSOURCE' 

//SYSIN  DD  * 

REAL  YB(900) ,XB(900) 

L0GICAL*1  SYMA/ ' / , SYMB/ ' , ' / , SYMC/ ' . ' / 

REAL*4  D(6 ,600) 

INTEGER* 2  A(600) ,B(600) ,C(600) 

LOGICAL* 1  LABEL(40, 2) 

***************************************************************** 

C  INFINITE  BEAM  ON  ELASTIC  FOUNDATION 

***************************************************************** 

C  SHEAR  AND  FLEXURAL  WAVES  BASED  ON  DYNAMIC  LAWS 
C  UNITS  IN  LB-IN-SEC  SYSTEM 

C  POSITIVE  SIGN  COWENTION  LISTED  AFTER  DEFINITION 
C 


c 

AIM 

-  ANGULAR  IMPULSE  DUE  TO  MOMENT 

LB-IN-SEC 

c 

AIN 

-  AREA  CONTRIBUTING  TO  DYNAMIC  INERTIA 

IN**2 

c 

AIV 

-  ANGULAR  IMPULSE  DUE  TO  SHEAR 

LB -IN-SEC 

c 

AS 

-  EFFECTIVE  SHEAR  CARRYING  AREA 

IN**2 

c 

Cl 

-  DILATATION  WAVE  VELOCITY 

IN/ SEC. 

c 

C2 

-  SHEAR  WAVE  VELOCITY 

IN/ SEC. 

c 

DEN 

-  WEIGHT  DENSITY  OF  BEAM  MATERIAL 

LB/CU  IN. 

c 

DT 

-  TIME  INCREMENT 

SEC. 

c 

DS 

-  DISTANCE  SCALE  FACTOR 

IN. 

c 

DTD 

*  TIME  INCREME NT -DI LA TATIONAL  WAVE 

SEC. 

c 

DTR 

-  TIME  INCREMENT-SHEAR  WAVE 

SEC. 

c 

DX 

-  CELL  LENGTH 

IN. 

c 

EM 

-  YOUNG'S  MODULUS  OF  ELASTICITY 

PSI. 

c 

ESR 

-  FOUNDATION  SPRING  CONSTANT 

PSI. 

c 

FMT 

-  PRINT  FORMAT  FOR  X-T  DIAGRAMS 

c 

G 

-  SHEAR  MODULUS  OF  RIGIDITY 

PSI. 

c 

IDE  NT 

-  TITLE  OF  RUN 

c 

IM 

-  NO.  OF  CELLS  INTO  WHICH  BEAM  IS  DIVIDED 

c 

INI 

-  MOMENT  OF  INERTIA  WHICH  CONTRIBUTES  RESISTANCE 

c 

TO  BENDING 

IN**4 

c 

INB 

-  MOMENT  OF  INERTIA  WHICH  CONTRIBUTES  RESISTANCE 

c 

TO  DYNAMIC  INERTIA 

IN*  *4 

c 

IPP 

-  0  PLOT  THE  DATA 

c 

IPP 

-  1  PRINT  THE  DATA 

c 

IPP 

-  2  PRINT  AND  PLOT  THE  DATA 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


IPR  -  NO.  OF  COLUMNS  OF  PRINTOUT  IN  A  PAGF 
KKM  »  TIME  CONTROL  INDEX  FOR  THF.  FIRST  RESPONSF.-DIST.  CURVE 

KKM1  =*  TIME  CONTROL  INDEX  FOR  THE  SECOND  RF.SPONSE-DIST.  CURVE 

KKM2  -  TIME  CONTROL  INDEX  FOR  THE  THIRD  RESPONSE-DIST.  CURVE 

KM  *  TIME  CONTROL  INDEX  FOR  TERMINATING  COMPUTATION 


KS  =  SHEAR  CORRECTION  FACTOR 

L  =*  LENGTH  OF  BEAM  IN. 

LIQ  -  LINEAR  IMPULSE  DUE  TO  FOUNDATION  REACTION  LB-SEC. 

LIQQ  -  LINEAR  IMPULSE  DUE  TO  DISTRIBUTED  LOAD  LB-SEC. 

LIV  -  LINEAR  IMPULSE  DUE  TO  SHEAR  LB-SEC. 

MM  -  SCALED  VALUE  OF  MOMENT 

MOM  -  INTERNAL  BENDING  MOMENT  -  CW  ON  LEFT  IN-LB. 

MS  »  MOMENT  SCALE  FACTOR  IN-LB. 

MO  -  GIVEN  APPLIED  MOMENT  AT  X-0  IN- LB. 

NAME  -  NAME  OF  MATERIAL  OF  BODY 

NU  »  POISSON'S  RATIO 

OMEGA  -  VELOCITY  OF  ROTATION  OF  CROSS-SECTION, 

D(PSI) /DT  -  CLOCKWISE  RAD/SEC. 

OMGS  =  ANGULAR  VELOCITY  SCALE  FACTOR  RAD/SEC. 

OMOM  =  SCALED  VALUE  OF  ANGULAR  VELOCITY 
PS I  -  ROTATION  OF  CROSS-SECTION  ABOUT  THE  NEUTRAL 

AXIS  -  CW  ON  LEFT  RAD. 

PS IPS  -  SCALED  VALUE  OF  ROTATION 

PSIS  -  ROTATION  SCALE  FACTOR  RAD. 

Q  -  INTENSITY  OF  UNIFORMLY  DISTRIBUTED  EXTERNAL 

LOADING  -  DOWNWARD  LB/ IN. 

QQO  -  INPUT  VALUE  OF  DISTRIBUTED  LOAD  LB/IN. 

RHO  -  MASS  DENSITY  OF  BEAM  MATERIAL  LB-SEC**2/IN** 

T  -  TIME  SEC. 

TAU1  -  DAMPING  TIME  CONSTANT  FOR  OMEGA  SEC. 

TAU2  -  DAMPING  TIME  CONSTANT  FOR  VELOCITY  SEC. 

TD  »  TIME  CLOCK  WHICH  REGULATES  DILATATIONAL  WAVE  SEC. 

TL  -  THE  PREVIOUS  TD  OR  TR  SEC. 

TLB  -  ONE  TIME  SEGMENT  BEFORE  TL  SEC. 

TR  -  TIME  CLOCK  WHICH  REGULATES  SHEAR  WAVE  SEC. 

TS  -  TIME  SCALE  FACTOR  SEC. 

TT  -  SCALED  VALUE  OF  TIME 

T1  -  DURATION  OF  MOMENT  PULSE  AT  X-0  SEC. 

T3  -  DURATION  OF  VELOCITY  PULSE  AT  X-0  SEC. 

V  -  VERTICAL  SHEAR  ON  A  CROSS-SECTION  -UP  ON  LEFT  LB. 

VEL  -  VELOCITY  OF  DEFLECTION,  D(Y)/DT  -  DOWNWARD  IN/SF.C. 

VELS  -  LINEAR  VELOCITY  SCALE  FACTOR  IN/SEC. 

VELVE  -  SCALED  VALUE  OF  LINEAR  VELOCITY 

VS  -  SHEAR  SCALE  FACTOR  LB. 

W  -  SCALED  VALUE  OF  SHEAR 

VO  -  GIVEN  APPLIED  SHEAR  AT  X-0  LB. 

X  -  AXIAL  COORDINATE  ALONG  LENGTH  OF  BEAM  MEASURED 

FROM  LEFT  END  IN. 

Y  -  DEFLECTION  -  DOWNWARD  IN. 

YS  -  DEFLECTION  SCALE  FACTOR  IN. 

YY  -  SCALED  VALUE  OF  DEFLECTION 


MAIN  PROGRAM 

DIMENSI0NX(400)  ,V  (400)  ,VEL(400)  ,Y(400)  ,PSI(400)  ,0MF.GA(400)  ,VV  (400) 


o  o 
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l ,7ELVE(400) , 0M0M(400 ) , YY(400) , PS I PS (400 ) , Q(400) , YB1 (900) , YB2(°00) 
REA LI MB , INI , L,M0M(400) , MO, MS , MM( 400) , NU ,KS , MASS , INERT 
L0CICAL*1IDENT(30) ,FMT(80) , NAME ( l 0) 

COMMONT,X,G,EM, L,MOM,V ,01 , C2 , Q,VEL, Y, RHO, DEN, PSI , OMECA, INB , INI ,AIN 
1  ,  AS ,  KS  , MO  ,7 0 ,  T1  ,  T3  ,  DT ,  DX ,  TS  ,  MS  ,V  S  ,7 ELS  ,  OMGS  ,  TT ,V'»  ,  MM  ,VELV E  ,  OMOM ,  TM 
l , KM , IPR, MASS , INERT, KKM , PS  IPS , PS  IS , YS , IMO , TAU1 , TAU2 , 000 , ESR , BETA 
1 , KKM1 , KKM2 
1  READ  801 , IDE NT, FMT 

READ  802,NAME,DEN,EM,L,NU,DX,AIN,KS 
IF(DEN)  500,500,3 

3  READ  803  ,  MO  ,V0 ,  T1 ,  T3 ,  TS  ,V  S  ,VELS  ,  YS 

READ  804  ,  KM  ,  IPR ,  INB ,  MS ,  PS  IS  ,  OMGS  ,  TAU1 ,  TAU2 
READ  805 , ESR, KXM, KKM1 , KKM2 ,  IPP , jS 
G-F.M/ ( 2 . 0*(  1 .  O+NU)  ) 

RH0-DEN/386.0 

AS«KS*AIN 

INI-INB 

C1»(((SM*INB)/(RH0*INI))**0.5) 

C2-( ( (AS*G)/(RHO*AIN) )**0. 5) 

MASS«RHO*AIN*DX 

INERT«RHO*INI*DX 

BETA- (ESR/ (4 .0*EM*INB) )**0. 25 

DT-DX/ (Cl ) 

GAMMA-EM/ ( G*KS ) 

C  DIVISION  OF  BEAM  INTO  ELEMENTS 
X(l)-0.0 
DO  101  1-1,400 
X(I+1)-X(I)+DX 

IF(ABS(X(I+l)-L)-DX/2.0)  10,10,101 

10  IM-I 

GO  TO  11 
101  CONTINUE 

1 1  IMO-IM+1 

PRINT  900 , IDENT 

PRINT  901, NAME, DEN, RHO, EM, G, ESR, NU 
PRINT  902 , L, AS , AIN , INB , INI , TAU1 , TAU2 
PRINT  903 ,C1 ,C2 
PRINT  904 ,M0,V0, T1 , T3 

PRINT  905 , DT, DX, TS  ,VS  ,VELS ,OMGS ,MS, PSIS , YS 
PRINT  906, IM, KM, IPR 

C  PRINTOUT  OF  CELL  BOUNDARY  COORDINATES 
PRINT  9266 

PRINT  907 , ( I , X(I ) , 1-1 , IMO) 

DO  909  M-11,13 

909  REWIND  M 

DO  910  M-20,22 

910  REWIND  M 
CALL  BEAM 
*********** 

PRINT  INSTRUCTIONS 
IF(IPP.EO.O)  GO  TO  500 
IA-1 
IM2-IPR 

30  DO  912  M-11,13 
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912  REWIND  M 

DO  913  M-20,22 

913  REWIND  M 
IF(TM+1-IM2)  24.24,25 

24  IM2-IM+1 

C  X-T  DIACRAM  FOR  MOMENT 
25  PRINT  908 
DO  104  K-l ,KM 

READ  (11)  TT,MM(1),VV(1),(MM(I),I-1 ,IM0) 

104  PRINT  FMT.TT,MM(1),VV(1),(MM(I),I-IA,IM2) 

C  X-T  DIAGRAM  FOR  SHEAR 
PRINT  1909 
DO  105  K-l, KM 

READ  (12)  TT,MM( 1 )  ,VV  (1 ) ,  (VV  (I) , 1-1 , IMO) 

105  PRINT  FMT,TT,MM(1),VV(1),(VV(I),I-IA,IM2) 

C  X-T  DIAGRAM  FOR  LINEAR  VELOCITY 

PRINT  1910 
DO  106  K-l, KM 

READ  (20)  TT.MM(l)  ,VV  ( 1 )  ,  (V ELVE (I  )  ,  1-1 ,  IMO) 

106  PRINT  FMT,TT,MM(1)  ,W(1)  ,(VELVE(I)  ,I-IA,IM2) 

C  X-T  DIAGRAM  FOR  DEFLECTION 

PRINT  911 
DO  107  K-l, KM 

READ  (21)  TT,MM( 1 )  ,W  (1),(YY(I),I-1,  IMO) 
PRINT  FMT,  TT,MM( 1 ) ,W (1 ) ,(YY(I ) , I -I A, IM2 ) 

107  CONTINUE 

C  X-T  DIAGRAM  FOR  ANGULAR  VELOCITY 

PRINT  1912 
DO  108  K-l, KM 

READ  (22)  TT,MM( l ) ,W (1 ) ,(OMOM(I ) , 1-1 , IMO) 
108  PRINT  FMT, TT,MM( 1 )  ,W  (1 ) ,(OMOM( I ) , I-IA, IM2) 

C  X-T  DIAGRAM  FOR  ANGLE  OF  ROTATION 
PRINT  4913 
DO  4108  K-l, KM 

READ  (13)  TT,MM( 1 ) ,VV  (1 ) ,( PS IPS (I ) , I-l , IMO) 
4108  PRINT  FMT, TT,MM( l )  ,VV  (l ) , (PS IPS (I ) , I-IA, IM2) 
IF(IM-IM2)  500,500,37 
37  IA-IA+IPR 
IM2-IM2+IPR 
GO  TO  30 

IF(IPP.EQ.l)  GO  TO  501 

800  FORMAT(80A1) 

801  F0RMAT(80A1/80AI) 

802  F0RMAT(10AI ,F10.2,E10.2,3F10.2,E10.2,F10.2) 

803  F0RMAT(8E10. 2) 

804  F0RMAT(2I10, 6E10, 2) 

805  F0RMAT(E10. 2,4I10,F10.2) 

900  F0RMAT(1HD,80A1//) 


901  FORMAT( 2  6H0MATERIAL  -  ,10A1 

/ 

1  , 24H0WEIGHT  DENSITY  DEN  -  ,F8.3,20H  LBS/IN**3  / 

2  , 24H0MASS  DENSITY  RHO  -  ,F8.6,20H  LB/SEC**2/IN**4  / 

3  , 24H0Y0UNGS  MODULUS  E  »  ,E8.2,20H  PS I  / 

4  , 24HOMOD .  OF  RIGIDITY  G  -  ,E8.2,20H  PSI  / 


5  , 24H0M0D.  OF  SUB.  RF.AC.ESR  -  ,E8.2,20H  PSI  / 

6  '  , 24HOPOISSONS  RATIO  NU  =  .F8.2.20H  /) 

902  F0RMAT(24H0LENGTH  OF  BEAM  L  -  .F8.3.20H  IN  / 

1  24HOSHEAR  AREA  AS  =  ,F8.2,17H  IN**2  / 

2  24H0INERTIA  AREA  AIN  -  ,F8.2,17H  IN**2  / 

3  24H0BEND.  MOM.  IN.  INB  -  ,F8.2,17H  IN**4  / 

4  24H0INERTIA  MOM.  IN.  TNI  -  ,F8.2,17H  IN**4  / 

5  24H0  TAU1  TAU1  -  ,E8.2,17H  1/SEC  / 

4  24HO  TAU2  TAU2  -  ,E8.2,17H  l/SEC  /) 

903  FORMAT(24HOMOMENT  WAVE  Cl  -  ,F8.0,20H  IN/SEC  / 

1  24HOSHEAR  WAVE  C2  -  ,F8.0,20H  IN/SEC  /) 

904  F0RMAT(24H0X-0  MOMENT  PULSE  MO  =  .F8.2.20H  IN-LB  / 

1  24H0X-0  VELOC  PULSE  VO  »  ,F8.2,20H  LB  / 

2  24H0X-0  MOM.  PULSE  DUR.T1*  ,E8.2,20H  SEC  / 

3  24H0X-O  VEL  PULSE  DUR.T3  -  ,E8.2,20H  SEC  /) 

905  FORMAT(24HOTIME  INCREMENT  DT  -  ,E8.2,20H  SEC  / 

1  24H0CELL  LENGTH  DX  -  ,F8.5,20H  IN  / 

2  24H0TIME  SCALE  FACTOR  TS  »  ,E8.2,20H  SEC  / 

3  24H0SHEAR  SCALE  FACTOR  VS-  ,E8.2,20H  LB  / 

4  24HWEL  SCALE  FACTOR  VELS-  ,E8.2,20H  IN/SEC  / 

5  24H0ANG  VEL  SCALE  FAC  OMGS-  ,E8.2,20H  RAD/SEC  / 

6  24H0M0MENT  SCALE  FACTOR  MS-  ,E8.2,20H  IN-LB  / 

7  24H0ANG  SCALE  FACTOR  PSIS-  ,E8.2,20H  RAD  / 

3  24H0DEF  SCALE  FACTOR  YS-  ,E8.2,20H  IN  /) 

906  F0RMAT( 24HONO .  OF  CELLS  IM  -  ,I4,17H  / 

1  24H0TOT.  TIME  INTERVALS  KM-  ,I4,17H  / 

2  24HONO.  COLS.  PRINTED  IPR-  ,14,  ////) 

907  F0RMAT(I4,F8.4) 


908  F0RMAT(1H1,2X, 'TIME', 5X, 'PULSE', 24X, 'MOMENT  PROPAGATION'/) 

911  F0RMAT(1H1,2X, 'TIME', 5X, 'PULSE', 21X, 'DEFLECTION  PROPAGATION'/) 

1909  F0RMAT(1H1,2X, 'TIME', 5X, 'PULSE', 24X, 'SHEAR  PROPAGATION'/) 

1910  F0RMAT(1H1, 2X, 'TIME', 5X, 'PULSE' ,16X, 'LINEAR  VELOCITY  PROPAGATION' 

l/) 

1912  F0RMAT(1H1, 2X, 'TIME', 5X, 'PULSE' ,15X, 'ANGULAR  VELOCITY  PROPAGATION 
1'/) 

4913  F0RMAT(1H 1 , 2X , ' TIME' , 5X , 'PULSE ' , 2 3X , ' ROTATION  PROPAGATION' / ) 

9266  F0RMAT(1HS ,2X,'I' ,3X,'X(I)' ,/) 

500  REWIND  21 

DO  127  K-l.KM 

READ  (21)  TT,MM(1) ,VV  (1) ,(YY(I) ,1-1 , IMO) 

XB(K)-TT 
YB(K)-YY(1 ) 

YB1(K)-YY(4) 

YB2(K)-YY(7) 

127  CONTINUE 

CALL  INITQ(A, B,C,D, 600) 

CALL  STSWQ(4662 , 30) 

READ(5 , 800)  (LABEL(J , l ) , J-l , 40) 

READ ( 5 , 800 )  (LABEL(J , 2) , J-l , 40) 

C  SET  GRPHIT  ARGUMENT  VALUES. 

N-KM 
LOGX-O 
LOGY-O 
XAXIS-8 . 0 
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YAXIS-5.0 

XMIN-0.0 

NXDEC-5 

XINC-20.0 

YMIN— 10.0 

NYDEC-5 

YINC -4.0 

XSMIN-0 . 0 

XSMAX-100.0 

YSMIN— 10.0 

YSMAX-10.0 

NDECX-1 

NDECY-1 

HT-0. 15 

HTS-0.15 

ICON-O 

CALL  PL0TIT(YB  ,  N,  LOGX ,  LOGY ,  XAXIS  ,  YAXIS  ,  XMIN,  XINC ,  NXDEC , 

1  YMIN, YINC , NYDEC , LABEL, XSMIN, XSMAX , YSMIN, YSMAX, XB , SYMA , ICON, HT , 

1  NDECX,NDECY,HTS) 

DO  130  I-1,N 
130  YB(I)-YB1(I) 

ICON-0 

CALL  SAMEP(YB, N,XB, SYMB, ICON) 

DO  71  I-l.N 
71  YB(I)-YB2(I) 

ICON-O 

CALL  SAMEP(YB , N, XB , SYMC , ICON) 

CALL  EPLOT 
REWIND  12 
DO  128  K-l.KM 

READ  (12)  TT,MM( 1 ) ,VV (1 ) ,(VV  (I ) , 1-1 , IMO) 

XB(K)-TT 

YB(K)-W(1) 

YB1(K)-W(4) 

YB2(K)-VV(7) 

128  CONTINUE 

CALL  INITQ(A,B,C,D, 600) 

CALL  STSWQ(4662 , 31 ) 

READ(5,800)  (LABEL(J, 1) ,J-1 ,40) 

READ(5,800)  (LABEL(J , 2 ) , J-l , 40) 

C  SET  GRPHIT  ARGUMENT  VALUES. 

N-KM 

LOGX-O 

LOGY-O 

XAXIS-8 . 0 

YAXIS-5.0 

XMIN-0.0 

NXDEC-5 

XINC-20.0 

YMIN— 10.0 

NYDEC-5 

YINC -4 . 0 

XSMIN-0. 0 

XSMAX-100.0 
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YSMIN— 10.0 

YSMAX- 10.0 

NDECX-1 

NDECY-1 

HT-0.15 

HTS-O. 15 

ICON-O 

CALL  PLOTIT(YB , N, LOGX , LOGY , XAXIS , YAXIS , XMIN , XI NC , NXDEC . 

1  YMIN, YINC , NYDEC , LABEL, XSMIN, XSMAX, YSMIN, YSMAX, XB , SYMA, ICON, HT, 
1  NDECX,NDECY,HTS) 

DO  72  I-1,N 

72  YB(I)-YB1(I) 

ICON-O 

CALL  SAMEP(YB , N,XB, SYMB, ICON) 

DO  73  1*1, N 

73  YB(I)-YB2(I) 

ICON-O 

CALL  SAMEP(YB, N,XB, SYMC, ICON) 

CALL  EPLOT 
REWIND  11 
DO  129  K-l.KM 

READ  (11)  TT,MM(1 )  ,VV  (1 ) , (MM(I ) ,1*1 , IMO) 

XB(K)*TT 
YB(K)-MM( 1 ) 

YB1 (K)-MM(4 ) 

YB2(K)-MM(7) 

129  CONTINUE 

CALL  INITQ( A, B , C , D , 600 ) 

CALL  STSWQ(4662 , 32 ) 

READ ( 5 , 800 )  (LABEL(J , 1 ) , J-l , AO) 

READ(5,800)  (LABEL(J , 2) , J-l , AO) 

C  SET  GRPHIT  ARGUMENT  VALUES. 

N-KM 

LOGX-O 

LOGY-O 

XAXIS-8 .0 

YAXIS-5 . 0 

XMIN-0.0 

NXDEC-5 

XINC-20. 0 

YMIN— 10.0 

NYDEC-5 

YINC-4.0 

XSMIN-0.0 

XSMAX-100.0 

YSMIN— 10. 0 

YSMAX-10.0 

NDECX-1 

NDECY-1 

HT-O. 15 

HTS-0.15 

ICON-O 

CALL  PLOTIT(YB , N, LOGX, LOGY , XAXIS , YAXIS , XMIN, XI NC , NXDEC, 

1  YMIN , YINC , NYDEC , LABEL, XSMIN , XSMAX , YSMIN , YSMAX , XB , SYMA , ICON, HT , 
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1  NDECX , NDECY , HTS ) 

DO  74  I-l.N 
74  YB(I)-YB1(I) 

ICON-O 

CALL  SAMEP(YB, N, XB, SYMB, ICON) 

DO  75  I-l.N 
75  YB(I )-YB2(I ) 

ICON-O 

CALL  SAMEP(YB,N,XB,SYMC, ICON) 

CALL  EPLOT 
REWIND  21 
DO  131  K-l.KM 

READ  (21)  TT,MM( 1 )  ,W  (1),(YY(I),I-1, IMO) 

IF(K. NE.KKM. AND.K. NE . KKM1 . AND. K. NE. KKM2)  GO  TO  131 

DO  231  I-l.IMO 

XB(I)-X(I) 

IF(K.EQ.KKM)  YB(1)-YY(I) 

IF(K.EQ.KKMl)  YB1(I)-YY(T) 

IF(K.EQ.KKM2)  YB2(I)-YY(I) 

231  CONTINUE 
131  CONTINUE 

CALL  INITQ(A,B,C,D,600) 

CALL  STSWQ(4662,33) 

READ(5 , 800)  (LABEL(J.l) ,J-1 ,40) 

READ(5 , 800)  (LABEL(J,2) ,J-1 ,40) 

C  SET  GRPHIT  ARGUMENT  VALUES. 

N-IMO 

LOGX-O 

LOGY-O 

XAXIS-8.0 

YAXIS-5 . 0 

XMIN-0.0 

NXDEC-5 

XINC-O. 2*DS 

YMIN— 10.0 

NYDEC-5 

YINC-4.0 

XSMIN-0 . 0 

XSMAX-DS 

YSMIN—  1 0.0 

YSMAX-10. 0 

NDECX-1 

NDECY-l 

HT-0.15 

HTS -0.15 

ICON-O 

CALL  PLOTIT(YB , N, LOGX , LOGY . XAXIS , YAXIS , XMIN, XINC , NXDEC , 

1  YMIN, YINC.NYDEC, LABEL, XSMIN.XSMAX, YSMIN, YSMAX, XB , SYMA, ICON, HT, 
1  NDECX, NDECY, HTS) 

DO  76  1-1 , N 
76  YB(I)-YBl(I) 

ICON-O 

CALL  SAMEP(YB,N,XB, SYMB, ICON) 

DO  77  I-l , N 
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77  YB(I)-YB2(I) 

ICON-0 

CALL  SAMEP(YB , N,XB , SYMC , ICON) 

CALL  EPLOT 
REWIND  12 
DO  132  K-l.KM 

READ  (12)  TT,MM(1) ,VV (1) ,(VV  (I) ,1-1 ,IMO) 

IF(K. NE.KKM. AND.K. NE . KKM1 . AND. K. NE . KKM2)  GO  TO  132 

DO  232  I-l.IMO 

XB(I)-X(I) 

IF(K.EQ.KKM)  YB(IHV(I) 

IF(K.EQ.KKMl)  YB1(I)-W(I) 

IF(K.EQ.KKH2)  YB2(I)-VV(I) 

232  CONTINUE 
132  CONTINUE 

CALL  INITQ( A , B , C , D , 600 ) 

CALL  STSWQ(4662 , 34) 

READ(5 , 800 )  (LABEL(J, 1) , J«1 ,40) 

READ (5 ,800)  (LABEL( J , 2) , J-l , 40) 

C  SET  GRPHIT  ARGUMENT  VALUES. 

N-IMO 

LOGX-O 

LOGY-O 

XAXIS-8 .0 

YAXIS-5. 0 

XMIN-0.0 

NXDEC-5 

XINC-0.2*DS 

YMIN— 10.0 

NYDEC-5 

YINC-4.0 

XSMIN-O. 0 

XSMAX-DS 

YSMIN—  1 0.0 

YSMAX-10. 0 

NDECX-1 

NDECY-1 

HT-0.15 

HTS-O. 15 

ICON-O 

CALL  PLOTIT(YB , N, LOGX, LOGY , XAXIS , YAXIS , XMIN, XINC , NXDEC, 
l  YMIN, YINC , NYDEC , LABEL, XSMIN , XSMAX , YSMIN , YSMAX , XB , SYMA , ICON , HT , 
1  NDECX , NDECY , HTS ) 

DO  78  1-1, N 
78  YB(I)-YBKI) 

ICON-O 

CALL  SAMEP (YB , N, XB , SYMB , ICON) 

DO  79  1-1 , N 
79  YB( I )-YB2(I ) 

ICON-O 

CALL  SAMEP(YB,N,XB, SYMC, ICON) 

CALL  EPLOT 
REWIND  11 
DO  133  X-l.KM 
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READ  (11)  TT,MM(1),VV(1),(MM(I),I-1,IM0) 

IF(K. NE.KKM.  AND.K. NE . KKM1 • AND. K. NE.KKM2)  GO  TO  133 

DO  233  I-l.IMO 

XB(I)-X(I) 

IF(K.EQ.KKM)  YB(I)-MM(I) 

IF(K.EQ.KKMl)  YBl(I)-MM(I) 

IF(K.EQ.KKM2)  YB2(I)-MM(I) 

233  CONTINUE 
133  CONTINUE 

CALL  INITQ(A, B,C,D,600) 

CALL  STSWQ(4662 , 35) 

READ ( 5 , 800 )  (LABEL(J, 1) , J-l ,40) 

READ(5 , 800 )  (LABEL(J , 2) , J-l ,40) 

C  SET  GRPHIT  ARGUMENT  VALUES. 

N-IMO 

LOGX-O 

LOGY-O 

XAXIS-8.0 

YAXIS-5.0 

XMIN-0.0 

NXDEC-5 

XINC-0.2*DS 

YMIN— 10.0 

NYDEC-5 

YINC-4.0 

XSMIN-0.0 

XSMAX-DS 

YSMIN— 10.0 

YSMAX-10.0 

NDECX-1 

NDECY-l 

HT-0.15 

HTS-0.15 

ICON-O 

CALL  PLOTIT(YB , N, LOCX , LOGY , XAXIS, YAXIS , XMIN, XI NC , NXDEC , 
l  YMIN, YINC , NYDEC , LABEL, XSMIN , XSMAX, YSMIN, YSMAX, XB , SYMA, ICON, HT , 
1  NDECX , NDECY , HTS ) 

DO  81  I-1,N 
81  YB(I)-YB1(I) 

ICON-O 

CALL  SAMEP(YB , N,XB, SYMB, ICON) 

DO  82  1-1 , N 
82  YB(I)-YB2(I) 

ICON-O 

CALL  SAMEP(YB,  N, XB.SYMC, ICON) 

CALL  EPLOT 
YB( 1 )-Q. 0 
YB1(1 )-0.0 
YB2(l )-0.0 
REWIND  21 
DO  136  K-l.KM 

READ  (21)  TT , MM( 1 ) ,VV (1),(YY(I),I-1, IMO ) 

IF(K.  NE.KKM. AND.K. NE.KKM1 .AND.K. NE .KKM2 )  GO  TO  136 
DO  236  I-l.IM 
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XB(I)-X(I) 

IF(K.EQ.KKM)  YB( 1+1 )-(YY(I+l )-YY(I) ) /DX 
IF( K .  EQ .  KKM1 )  YB  1  ( 1+1  )=*  (  YY  ( 1+1  )-YY  ( I )  )  /DX 
IF(K.EQ.KKM2)  YB2(I+1 )»(YY(I+1 )-YY(I))/DX 
CONTINUE 
CONTINUE 

CALL  INITQ(A,B,C,D,600) 

CALL  STSWQ(4662 , 36) 

READ ( 5 , 800 )  (LABEL(J , 1 ) , J-l , 40) 

READ( 5 , 800)  (LABEL(J,2) ,J“1 ,40) 

SET  GRPHIT  ARGUMENT  VALUES. 

N-IMO 

LOGX-O 

LOGY-O 

XAXIS-8.0 

YAXIS-5.0 

XMIN-0.0 

NXDEC-5 

XINC-0.2*DS 

YMIN— 1.0 

NYDEC-5 

YINC-0.4 

XSMIN-0.0 

XSMAX-DS 

YSMIN— 1.0 

YSMAX-1 . 0 

NDECX-1 

NDECY-1 

HT-0.15 

HTS-0.15 

ICON-O 

CALL  PLOTIT( YB , N, LOGX , LOGY , XAXIS , YAXIS , XMIN, XINC , NXDEC , 

1  YMIN , YI NC , NYDEC , LABEL, XSMIN, XSMAX .YSMIN , YSMAX , XB , S YMA , IC  ON,  HT , 

1  NDECX, NDECY.HTS) 

DO  84  I-l.N 
YB(I)-YB1(I) 

ICON-O 

CALL  SAMEP(YB, N,XB, SYMB, ICON) 

DO  85  I-l.N 
YB(I)-YB2(I) 

ICON-O 

CALL  SAMEP(YB , N, XB , SYMC , ICON) 

CALL  EPLOT 

STOP 

END 

SUBROUTINE  BEAM 

DIMENSI0NX(400)  ,V  (400)  ,VEL(400)  ,Y(400)  ,PSI(400)  ,0MEGA(400)  ,\n/(400) 
1 ,VELVE(400) ,OMOM( 400) ,YY(400) ,PSIPS(400) ,AIM(400) ,AIV (400) ,Q(400) 
REALINB , INI , L,M0M(400) ,M0, MS ,MM(400) , NU.KS, MASS, INERT, LIV  (400) ,LIQ 
1(400) 

L0GICAL*1IDENT(80),FMT(80),NAME(10) 

COMMONT.X, G, EM, L,MOM,V , Cl , C2 , Q,VEL, Y,RHO, DEN, PS I , OMEGA, INB , INI , AIN 
1 , AS , KS , MO ,V 0 , T1 , T3 , DT , DX , TS , MS ,VS,VELS,OMGS,TT,VV  ,MM,VEIVE , OMOM,  IM 
1 , KM , IPR , MASS , INERT, KKM , PS IPS , PS IS , YS , IMO , TAU1 , TAU2 , QQO , ESR , BETA 


I ,KKM1 ,KKM2 
PRINT  921 

C  INITIAL  STATE  OF  BODY 
DO  119  I-l, IMO 
AIM(I)-0.0 
AIV(I)-0.0 
LIV (I)-O.O 
LIQ(I)-0.0 
LIQQ(I)-0.0 
M0M(I)-0.0 
0MEGA(I)-0.0 
PSI(I)«0.0 
V(I)»0.0 
VEL(I)»0.0 
119  Y(I)-0.0 

QQ-QQO 
DTD-DX/(C1) 

DTR-DX/(C2) 

T-0.0 
TD-0.0 
TL-0 . 0 
TLB -0.0 
TR-0.0 

DO  101  K-l.KM 
TT-T/TS 

C  ***INPUT  AT  X-0,  MOMENT  BOUNDARY  CONDITION*** 
IF(T-Tl)  110,111,111 

110  MOM(IMO)-MO 
GO  TO  1492 

111  MOM(IMO)-O.0 
1492  CONTINUE 

C  ***INPUT  AT  X-0,  SHEAR  BOUNDARY  CONDITION*** 

IF(T-T3)  113,114,114 

113  V(1)^0 
GO  TO  115 

114  V(l)-0.0 

115  CONTINUE 

C  ***REACTION  OF  ELASTIC  FOUNDATION*** 

DO  116  I-l.IM 

Q(I)»- ESR*Y(I)/(1+ESR*Y(I)/100.0) 

116  CONTINUE 

C  SCALING  FOR  PRINTOUT 

IZZ-IMO 

DO  103  I-i.IZZ 
MM(I)-MOM(I)/MS 
OMTM( I )-OMEGA( I )/OMGS 
PSIPS(I)-PSI(I)/PSIS 
VELVE(I)*?EL(I)/VELS 
W(I)-V(I)/7S 
YY(I)-Y<I)/YS 
103  CONTINUE 

WRITE  (11)  TT,MM(1 )  ,VV (1 ) ,(MM( I) , 1-1 , IMO) 
WRITE  (12)  TT , MM( 1 ) ,W (1 ) , (VY (I ) , I-l , IMO ) 
WRITE  (13)  TT,MM( 1 ) ,W (1 ) , (PSIPS(I ) , I-l , IMO) 


WRITE  (20)  TT,MM( 1 ) ,VV (l ) , (V  E  WE(I ) , I-l , IMO) 
WRITE  (21)  TT , MM( 1 )  , VV  ( 1 ) , ( YY ( I ) , I -1 , IMO ) 

WRITE  (22)  TT,MM(  1 )  ,W  (1 )  , ( OMOM(I )  ,  I-l ,  IMO) 

12  IF( ABS(TD-(T+DT) )-0 . IE-10)  19,19,89 

89  IF(TD-TR)  15,15,16 

16  IF(TD-H)TD-TR-DTR)  15,15,13 
C  PROPAGATION  PROCEDURE 

C  THE  Cl  WAVE 

15  DO  108  1-2, IMO 

C  ANGULAR  IMPULSE  DUE  TO  MOMENT 

AIM( I )- (MOM( I-l )-MOM( I ) )*DTD+AIM( I ) 

C  ANGULAR  IMPULSE  DUE  TO  SHEAR 

AW  (I )-(V  (I-l  )+V  (I)  ) *DX* ( TD+DTD-TL) / 2  .  O+AW  (I ) 
C  CUMULATE  LINEAR  IMPULSE  DUE  TO  THIS  SHEAR 

LW  (I-l  )-LW  (I-l  )+(V  (I)-V  (I-l))*(TtH-DTD-TL) 

C  DAMPING,  FIRST  HALF 

IF(TAU1-1 .0E+40)  930,928,928 
930  DOMG«(-1.0)*OMEGA(I)*(DTD/TAU1) 

OMEGA( I )-OMEGA( I )+(D0MG*0 . 5 ) 

C  IMPULSE-MOMENTUM  ACROSS  Cl 
928  DOM-(ATM(I)+AW  (I)) /INERT 
OMEGA( I )-OMEGA( I  )+DOM 
C  BOUNDARY  CONDITION  FOR  OMEGA 
0MEGA(l)-0.0 
PSI(l)-0.0 

C  DAMPING,  SECOND  HALF 

IF(TAU1-1 . 0E+401  9113,815,815 
9113  DOMG»(-l.O)*OMEGA(I)*(DTD/TAUl) 

OMEGA( I )»OMEGA( I )+(DOMG*Q . 5 ) 

C  INITIALIZE  ANGULAR  IMPULSES  AFTER  USE 
815  AW  (I  )-0 . 0 
AIM( I )-0 . 0 
108  CONTINUE 

TD-TIH-DTD 
TLB-TL 
TL-TD 

DO  107  I-l.IM 

DALPH-( (OMEGA( 1+1 )-OMEGA( I ))*DTD)/DX 
C  CONSTITUTWE  EQUATION  FOR  INCREMENTAL  MOMENT 

DMOM-(-l . Q)*EM*INB*DALPH 
C  CUMULATE  MOMENT 

MOM(I)«MCM(I)+DMOM 

C  CUMULATE  SHEAR  (ENTIRE  OMEGA  PORTION) 

V (1+1 )-V  (1+1 )-AS*G*OMEGA( 1+1 )*DTD 
C  ROTATION 

PSI(I+1 )-PSI(I+l )+(OMEGA(I+l )*DTD) 

107  CONTINUE 

GO  TO  12 

C  THE  C2  WAVE 

13  DO  104  I-l.IM 

C  LINEAR  IMPULSE  DUE  TO  FOUNDATION  REACTION 

LIQ(I )-Q(I )*DX*DTR+LIQ(I ) 

C  LINEAR  IMPULSE  DUE  TO  SHEAR 

LW  ( I )- (V  ( 1+1  )-V  ( I )  ) * (  TR+OTR-TD+L W  ( I ) 
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C  CUMULATING  ANGULAR  IMPULSE  DUE  TO  SAME  SHEAR 

AIV  (1+1  )-AIV  (1+1  )+0.  5 *(V  (I  )-+V  (1+1 ))*DX*(TR+DTR-TL) 

C  DAMPING,  FIRST  HALF 

IF(TAU2-1 . 0E+40 )  925,931,931 
925  DVELG-(-l .0)*VEL(I)*(DTR/TAU2) 

VEL(I)-VEL(I)+(DVELG*0. 5) 

C  IMPULSE-MOMENTUM  ACROSS  C2 
931  DVEL-(LIV (I )+LIQ(I ) ) /MASS 
VEL(I)-VEL(I)+DVEL 

C  BOUNDARY  CONDITION  FOR  VELOCITY 
VEL(IM+1 )-0.0 
C  DAMPING,  SECOND  HALF 

IF(TAU2-l . 0E+40)  927,818,818 
927  DVELG-(-l.O)*VEL(I)*(DTR/TAU2) 

VEL(I  )-VEL(I  )+(DVELG*0 .5) 

C  INITIALIZING  LINEAR  IMPULSES  AFTER  USE 
818  LIV(I)-0.0 

LIQ(I)-0.0 
104  CONTINUE 
TR-TR+DTR 
TLB-TL 
TL-TR 

DO  102  I-l.IM 
C  LINEAR  STRAIN 

DSTRN-((VEL(I+1)-VEL(I))*DTR)/DX 

C  CONSTITUTIVE  EQUATION  FOR  INCREMENTAL  SHEAR  (VELOCITY  PORTION) 
DV-AS*G*DSTRN 
C  CUMULATE  SHEAR 
V(I+1)-V(I+1)+DV 
C  DEFLECTION 

Y(I)«Y(I)+(VEL(I)*DTR) 

102  CONTINUE 
GO  TO  12 
19  CONTINUE 
T-T+DT 

101  CONTINUE 
ENDFILE  11 
ENDFILE  12 
ENDFILE  13 
ENDFILE  20 
ENDFILE  21 
ENDFILE  22 
921  FORMAT(lHl) 

500  RETURN 
END 

//DATA.FT11F001  DD  UNIT-SYSDA , SPACE-(TRK, (20, 10) ,RLSE) 

//DATA. FT12F001  DD  UNIT-SYSDA, SPACE-(TRK, (20, 10) ,RLSE) 

/ /DATA . FTl 3F001  DD  UNTT»SYSDA, SPACE* (TRK, (20,10) ,RLSE) 

/ /DATA . FT20F001  DD  UNIT-SYSDA, SPACE-(TRK, (20, 10) ,RLSE) 
//DATA.FT21F001  DD  UNIT-SYSDA, SPACE-(TRK, (20, 10) ,RLSE) 
//DATA.FT22F001  DD  UNIT-SYSDA, SPACE-(TRK, (20, 10) ,RLSE) 

//DATA. INPUT  DD  * 

BOLEY  CHECK 

(1HS ,F6. 2 , 1X,F5. 1 , 1X,F5. 1 , IX, 22F5. 2) 


CONCRETE 

0.0868 

30.00E+5 

48.00 

0.20 

1.00 
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4.0E4O 

0.833 

O.OOE+O 

-100.0E+0 

3 .  OOF.  40 

3  .  OOF. 40 

1.00E-4 

1 2 . OOE+O 

1  .OOE-2 

1.50E-4 

900 

22 

5.33E40 

50. 00E40 

1 .00F.-5 

1  .OOF. -2 

3.00E-5 

3.00E-5 

1.00E+4 

50 

100 

150 

50.00 

TIME  IN  1 . 0*E-4  SEC 
BEAM  DEFLECTION  IN  1.5*F.-4  INCHES 
TIME  IN  1 . 0*E-4  SEC 
BEAM  SHEAR  IN  12.0*LB. 

TIME  IN  1 . 0*E-4  SEC 
BEAM  MOMENT  IN  50.0*LB-IN. 

DISTANCE  X  IN  INCHES 

BEAM  DEFLECTION  IN  1.5*E-4  INCHES 

DISTANCE  X  IN  INCHES 

BEAM  SHEAR  IN  12.0*LB 

DISTANCE  X  IN  INCHES 

BEAM  MOMENT  IN  50.0*LB-IN. 

DISTANCE  X  IN  INCHES 
BEAM  SLOPE  IN  1.5*E-4 

/ /DATA . FT30F001  DD  DSN-MEN. PI 1 850 . A3B . PLOT, 

//  VOL-REF-MEN. P11850. A3B. LIB, DISP-(NEW, KEEP), 

//  SPACE-(TRK,(20,2),RLSE), 

/ /  DCB-( RECFM-FB , LRECL-80 , BLKSIZE-1 2960 ) 

//DATA. FT31F001  DD  DSN-MEN. PI 1850. A3B.PL0T1 , 

/  /  V  OL-REF-MEN. PI 1 850 . A3B . LIB , DISP- ( NEW , KEEP ) , 

//  SPACE- (TRK, (20,2) ,RLSE) , 

//  DCB-(RECFM-FB,LRECL-80, BLKSIZE-1 2960) 

/ /DATA. FT32F001  DD  DSN-MEN. PI 1850. A3B . PL0T2 , 

//  V OL-REF-MEN. PI 1850. A3B . LIB, DISP- (NEW, KEEP) , 

//  SPACE-(TRK,(20,2) ,RLSE) , 

//  DCB-(RECFM-FB,LRECL-80, BLKSIZE-1 2960) 

//DATA. FT33F001  DD  DSN-MEN. PI 1850. A3B . PL0T3 , 

//  VOL-REF-MEN. PI 1850. A3B. LIB, DISP-(NEW, KEEP) , 

//  SPACE-(TRK, (20,2) ,RLSE) , 

/ /  DCB-( RECFM-FB , LRECL-80 , BLKSIZE-1 2960 ) 

/ /DATA. FT34F001  DD  DSN-MEN. PI  1850. A3B . PL0T4 , 

/ /  V  OL-REF-MEN. PI 1 8  50 ■ A3B . LIB , DISP- (NEW , KEEP) , 

//  SPACE-(TRK,  (20,2) , RLSE) , 

/ /  DCB»( RECFM-FB , LRECL-80 , HLKS IZE-1 2960 ) 

/ /DATA. FT35F001  DD  DSN-MEN. P11850. A3B .PL0T5, 

//  VOL-REF-MEN.  PU850.  A3B. LIB, DISP-(NEW, KEEP)  , 

//  SPACE- (TRK, (20,2) ,RLSE) , 

/ /  DCB»( RECFM-FB , LRECL-80 , BLKSIZE-1 2960) 

/ /DATA. FT36F001  DD  DSN-MEN. P11850. A3B.PLOT6, 

//  VOL-REF-MEN. P11850.A3B. LIB, DISP-(NEW, KEEP) , 

//  SPACE- (TRK, (20,2) ,RLSE) , 

//  DCB-( RECFM-FB, LRECL-80, BLKSIZE-1 2960) 
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APPENDIX  B 

Computer  Program  for  Infinite  Plate 
on  Elastic  Foundation 


MICROCOPY  resolution  test  chart 

NSliONAl.  BURl‘1'  T»<  f-lNDARDS  I  “tv-  4 


i 
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//  EXEC  PCM-IEFBR14 

/ /D  DD  DSN-MEN. PU850.A3B. PLOT, DISP«(OLD, DELETE), 

//  VOL-REF-MEN. PI  1850. A3B . LIB 

//D  DD  DSN-MEN. PI 1850 . A3B . PL0T1 ,DISP-( OLD, DELETE) , 

//  VOL-REF-MEN. P11850.A3B. LIB 

// D  DD  DSN-MEN. PI  1850. A3B . PL0T2 , DISP-/OLD, DELETE), 

//  VOL-REF-MEN. PI 1850. A3B. LIB 

//D  DD  DSN-MEN. P11850.A3B.PLOT3,DISP»( OLD, DELETE) , 

//  VOL-REF-MEN. PI 1850. A3B . LIB 

// D  DD  DSN-MEN. PU850.A3B.PL0T4,DISP-(0LD, DELETE)  , 

//  V  OL-REF-MEN. P11850 . A3B . LIB 

//D  DD  DSN-MEN. PI  1850. A3B . PL0T5, DISP-(OLD, DELETE ) , 

//  V  OL-REF-MEN. PI 1 850 . A3B . LIB 

//D  DD  DSN-MEN. P11850. A3B . PL0T6, DISP-( OLD .DELETE) , 

//  VOL-REF-MEN. P11850.A3B. LIB 
//  EXEC  F GC LG, PARM-' NOSOURCE' 

//SYSIN  DD  * 

**************************************************************** 

C  CIRCULAR  PLATE  ON  ELASTIC  FOUNDATION 

**************************************************************** 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SHEAR  AND  FLEXURAL  WAVES  BASED  ON  DYNAMIC  LAWS 
UNITS  IN  LB-IN-SEC  SYSTEM 


AIMR 

ANGULAR  IMPULSE  DUE  TO  RADIAL  MOMENT 

LB -IN-SEC 

AIMTH 

ANGULAR  IMPULSE  DUE  TO  TANGENTIAL  MOMENT 

LB-IN-SEC 

AIQR 

* 

ANGULAR  IMPULSE  DUE  TO  TRANSVERSE  SHEAR 

LB-IN-SEC 

ALPHA 

m 

ANCULAR  STRAIN  -  D(PHI)/DR 

RAD/IN. 

CP 

« 

PLATE  VELOCITY 

IN/ SEC. 

C2 

m 

SHEAR  WAVE  VELOCITY 

IN/ SEC. 

D 

m 

FLEXXURAL  RIGIDITY 

LB-IN. 

DALPH 

m 

INCREMENTAL  ANGULAR  STRAIN 

RAD/IN. 

DEN 

m 

WEIGHT  DENSITY  OF  THE  PLATE  MATERIAL 

LB/CU  IN. 

DR 

m 

CELL  LENGTH 

IN. 

DSTRN 

a 

INCREMENTAL  LINEAR  STRAIN 

DS 

m 

DISTANCE  SCALE  FACTOR 

IN. 

DT 

m 

TIME  INCREMENT 

SEC. 

DTD 

m 

TIME  INCREME NT -DILATATIONAL  WAVE 

SEC. 

DTR 

m 

TIME  INCREMENT-SHEAR  WAVE 

SEC. 

EM 

m 

YOUNG'S  MODULUS  OF  ELASTICITY 

PSI. 

ESR 

m 

FOUNDATION  SPRING  CONSTANT 

PSI. 

FMT 

m 

PRINT  FORMAT  FOR  X-T  DIAGRAMS 

G 

■m 

SHEAR  MODULUS  OF  RIGIDITY 

PSI. 

H 

m 

THICKNESS  OF  THE  PLATE 

IN. 

IDE  NT 

m 

TITLE  OF  RUN 

IM 

m 

NO.  OF  CELLS  INTO  WHICH  PLATE  IS  DIVIDED 

IPP 

m 

0  PLOT  THE  DATA 

IPP 

m 

1  PRINT  THE  DATA 

IPP 

m 

2  PRINT  AND  PLOT  THE  DATA 

IPR 

M 

NO.  OF  COLUMNS  OF  PRINTOUT  IN  A  PAGE 

KKM 

m 

TIME  CONTROL  INDEX  FOR  THE  FIRST  RESPONSE¬ 

DIST.  CURVE 

KKM1 

m 

TIME  CONTROL  INDEX  FOR  THE  SECOND  RESPONSE 

-DIST.  CURVE 

KKM  2 

m 

TIME  CONTROL  INDEX  FOR  THE  THIRD  RESPONSE- 

DIST.  CURVE 

KM 

m 

TIME  CONTROL  INDEX  FOR  TERMINATING  COMPUTATION 

K2 

m 

SHEAR  CORRECTION  FACTOR 

o  o  a 


RAD/SEC. 
RAD /SEC. 


LB/IN. 

LB/IN. 

IN. 

LB-SEC**2/IN**4 

IN. 

IN. 


LIQR  -  LINEAR  IMPULSE  DUE  TO  TRANSVERSE  SHEAR  LB-SEC. 

MMR  -  SCALED  VALUE  OF  RADIAL  MOMENT 

MMRO  -  SCALED  VALUE  OF  R-RO  RADIAL  MOMENT  INPUT  IN- LB. 

MMTH  -  SCALED  VALUE  OF  TANGENTIAL  MOMENT  IN-LB. 

MR  -  RADIAL  BENDING  MOMENT  PER  UNIT  LENGTH  LB-IN/IN. 

MRS  -  RADIAL  MOMENT  SCALE  FACTOR  LB-IN/IN. 

MRO  -  GIVEN  APPLIED  RADIAL  MOMENT  AT  R-0  IN-LB. 

MTH  -  TANGENTIAL  BENDING  MOMENT  PER  UNIT  LENGTH  LB-IN/IN. 

MTHS  -  TANGENTIAL  MOMENT  SCALE  FACTOR  LB-IN/IN. 

NAME  -  NAME  OF  MATERIAL  OF  BODY 

NU  -  POISSON'S  RATIO 

OMEGA  -  VELOCITY  OF  ROTATION  OF  CROSS-SECTION, 

D(PSI)/DT  -  CLOCKWISE  RAD/SEC. 

OMGS  -  ANGULAR  VELOCITY  SCALE  FACTOR  RAD/SEC. 

OMOM  -  SCALED  VALUE  OF  ANGULAR  VELOCITY 

PHI  -  ROTATION  OF  THE  CROSS-SECTION  ABOUT  THE 

TANGENTIAL  AXIS  RAD. 

QQ  -  SCALED  VALUE  OF  SHEAR 

QQO  -  SCALED  VALUE  OF  R-RO  SHEAR  INPUT 

QR  -  TRANSVERSE  SHEAR  PER  UNIT  LENGTH  LB/IN. 

QRS  -  SHEAR  SCALE  FACTOR 

QRO  -  GIVEN  APPLIED  TRANSVERSE  SHEAR  AT  R-RO  LB/IN. 

R  -  RADIAL  COORDINATE  FROM  THE  CENTER  OF  PL  IN. 

RHO  -  MASS  DENSITY  OF  PLATE  MATERIAL  LB-SEC* *2 /II 

RL  -  OUTER  RADIUS  OF  THE  PLATE  IN. 

RO  -  INNER  RADIUS  OF  THE  PLATE  IN. 

STRAIN  -  LINEAR  STRAIN  -  0(W)/DR 

T  -  TIME  SEC. 

TAU1  -  DAMPING  TIME  CONSTANT  FOR  OMEGA  SEC. 

TAU2  -  DAMPING  TIME  CONSTANT  FOR  VELOCITY  SEC. 

TD  -  TIME  CLOCK  WHICH  REGULATES  DILATATIONAL  WAVE  SEC. 

TL  -  THE  PREVIOUS  TD  OR  TR  SEC. 

TR  -  TIME  CLOCK  WHICH  REGULATES  SHEAR  WAVE  SEC. 

TS  -  TIME  SCALE  FACTOR  SEC. 

TT  -  SCALED  VALUE  OF  TIME 

T1  -  DURATION  OF  MOMENT  PULSE  AT  R-0  SEC. 

T3  -  DURATION  OF  VELOCITY  PULSE  AT  R-0  SEC. 

VEL  -VELOCITY  OF  DEFLECTION,  D(W)/DT  -  DOWNWARD  IN/SEC. 

VELS  -  LINEAR  VELOCITY  SCALE  FACTOR  IN/SEC. 

VEUVE  -  SCALED  VALUE  OF  LINEAR  VELOCITY 

W  -  TRANSVERSE  DISP.  OF  THE  MIDPLANE  IN. 

WS  -  DEFLECTION  SCALE  FACTOR 

WW  -  SCALED  VALUE  OF  DEFLECTION 

Z  -  K2*C2  IN/ SEC. 

REAL  YB(900) ,XB(900) ,YB1(900) ,YB2(900) 

LOGICAL*l  SYMA/'"'/,SYMB/' ,'/ ,SYMC/' .'/ 

REAL**  D ( 6 , 600 ) 

INTEGER*2  A(600) ,B(600) ,C(600) 

L0GICAL*1  LABEL(40, 2) 

MAIN  PROCRAM 

DIMENSION  R(250),QR(250),VEL(250),W(250),0MEGA(25O),QQ(250), 
IVELVE(250)  ,0M0M(250)  ,WW(250)  ,U(250)  ,QQS(250) 

REAL  MR(250),MTH(250),NU, MRO, MRS, MTHS, MMR<250),MMTH(250), 
1MMR0.K2 


SEC. 

SEC. 

IN/ SEC. 
IN/ SEC. 


IN/ SEC. 
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L0GICAL*1  IDENT(80) ,FMT(80)  ,  NAME(IO) 

COMMON  T , R , G ,  EM , IMO , H , DI , RO , RL, MR , MTH , QR , CP , C2 , K2  ,V E L , W , RHO , 
1DF.N,  NU ,  OMEGA, MRO ,  QRO ,  T1 , T3  ,  DSTRN, DALPH , DT,  DR,  TS , MRS , 

2MTHS , QRS ,  V  E LS , OMGS , WS , TT , QQ , MMR , MMTH , V E LV E , OMOM , IM , KM , I PR , 
3WV,MMRO,QQO,QQS,Z,U,TAUl,TAU2,ESR,NAME,FMT,IDENT 
READ  801 , IDENT, FMT 
READ  802 , NAME, DEN, EM, RO, RL, NU, DR, H 
READ  803 , MRO , QRO , T1 , T3 , TAU1 , TAU2 
READ  803, TS, QRS ,VELS,WS, MRS, MTHS , OMGS 
READ  804, KM, IPR,KKM,KKM1,KKM2 
READ  805,ESR,DS , IPP 

800  F0RMAT(80A1 ) 

801  F0RMAT(80A1/80A1) 

802  FORMAT(10A1,F10.6,E10.2,3F10.2,F10.3,E10.2) 

803  F0RMAT(7E10. 4) 

804  F0RMAT(5I 10) 

805  F0RMAT(2E10. 2, 110) 

GEE-386.0 

C-EM/(2.0*(1.0+NU)) 

DI-(EM*(H**3  ))/((12.0)*(l  . 0-(NU**2  ))) 

RHO-DEN/GEE 

CP-(EM/ ( RHO*( 1 . 0-( NU**2 ) ) ) ) **0 . 5 

C2-(G/RH0)**0 . 5 

K2-( (0 . 76+(0 . 3*NU))**0 . 5) 

Z«K2*C2 

DT-DR/CP 

C  DISCRETATION  OF  PLATE  INTO  ELEMENTS 
R(1)-R0 

DO  101  1-1,250 
R(I+1 )«R(I)+DR 

IF(ABS(R(I+l )-RL)-DR/2 .0)  10, 10, 101 

10  IM-I 

GO  TO  11 
101  CONTINUE 

11  IMO-IM+1 
PRINT  900, IDENT 

900  F0RMAT(1HD,8QA1,/ /) 

PRINT  901 , NAME, DEN, RHO, EM, G, NU,DI 
PRINT  902 ,RO,RL,H 
PRINT  903,CP,C2,K2,Z 
PRINT  904,MR0,QR0,T1,T3 

PRINT  905 ,DT, DR, TS ,QRS ,VELS , OMGS .MRS , MTHS ,WS 
PRINT  906,IM,KM, IPR 
PRINT  9266 

9266  F0RMAT(1HS,2X,'I',3X,'R(I)',/) 

PRINT  907,(I,R(I) ,1-1 , IMO) 

MMRO-MRO/MRS 
QQO-QRO/QRS 
REWIND  11 
REWIND  12 
REWIND  13 
REWIND  20 
REWIND  21 
REWIND  22 
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CALL  PLATE 

IF(IPP.EQ.l)  GO  TO  500 
C  PRINT  INSTRUCTIONS 
IA-1 
IM2-IPR 

30  REWIND  II 
REWIND  12 
REWIND  13 
REWIND  20 
REWIND  21 
REWIND  22 

IF(IM+1-IM2)  24,24,25 

24  IM2-IM+1 

C  R-T  DIAGRAM  FOR  RADIAL  MEMENT 

25  PRINT  908 

DO  104  K-l.KM 

READ(ll)  TT.MMRO, QQO , (MMR(I ) , 1-1 , IMO) 
PRINT  FMT, TT.MMRO, QQO, (MMR(I) ,I-IA, IM2) 

104  CONTINUE 

C  R-T  DIAGRAM  FOR  TANGENTIAL  MOMENT 

PRINT  4913 
DO  4108  K-l.KM 

READ( 1 3 )  TT , MMRO , QQO , ( MMTH( I ) , 1-1 , IMO ) 
PRINT  FMT, TT.MMRO, QQO, (MMTR(I ) , I-IA, IM2) 
4108  CONTINUE 

C  R-T  DIAGRAM  FOR  SHEAR 

PRINT  1909 
DO  105  K-l.KM 

READ (12)  TT.MMRO, QQO, (QQ(I),I»1, IMO) 

PRINT  FMT, TT.MMRO, QQO, (QQ(I), I-IA, IM2) 

105  CONTINUE 

C  R-T  DIAGRAM  FOR  LINEAR  VELOCITY 

PRINT  1910 
DO  106  K-l.KM 

READ(20)  TT.MMRO, QQO, (VELVE(I) ,1-1 , IMO) 
PRINT  FMT, TT.MMRO, QQO, (VEIN E(I ) , I-IA, IM2) 

106  CONTINUE 

C  R-T  DIAGRAM  FOR  DEFLECTION 

PRINT  911 
DO  107  K-l.KM 

READ(21)  TT.MMRO, QQO, (WW(I),I-1, IMO) 

PRINT  FMT, TT.MMRO, QQO, (WW(I), I-IA, IM2) 

107  CONTINUE 

C  R-T  DIAGRAM  FOR  ANGULAR  VELOCITY 

PRINT  1912 
DO  108  K-l ,KM 

READ(22)  TT.MMRO, QQO, (OMOM(I).I-l, IMO) 
PRINT  FMT, TT.MMRO, QQO, (OMOM( I), I-IA, IM2) 

108  CONTINUE 

IF(IM-IM2)  9872,9872,37 
37  IA-IA+IPR 

IM2-IM2+IPR 
GO  TO  30 
9872  PRINT  913 
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C 


901 

F0RMAT( 2  6H0MATERIAL 

,  10A1 

/ 

1 

26H  WEIGHT  DENSITY 

DEN 

,  F8. 3 , 16H 

LBS/IN**3 

/ 

2 

26H  MASS  DENSITY 

RHO 

,  F8. 6, 16H 

LB/IPS**2/IN**4  / 

3 

26H  YOUNGS  MODULUS 

E 

,E8.2, 16H 

PSI 

/ 

4 

26H  MOD.  OF  RIGIDITY 

G 

,E8. 2, 16H 

PS  I 

/ 

5 

26H  POISSONS  RATIO 

NU 

,F8.2, 16H 

/ 

6 

26H  FLEXURAL  RIGIDITY 

D 

,E8. 2, 16H 

LB-IN 

/) 

902 

FORMAT(26HOINNER  RADIUS 

RO 

, F8. 4 , 16H 

IN 

/ 

1 

26H00UTER  RADIUS 

RL 

,F8.4, 16H 

IN 

/ 

2 

26H0PLATE  THICKNESS 

H 

,F8. 4, 16H 

IN 

/) 

903 

FORMAT( 2 6H0PLATE  VELOCITY 

CP 

,F10. 2, 16H 

IN/ SEC 

/ 

1 

26H  SHEAR  VELOCITY 

C2 

, F10. 2 , 16H 

IN/SEC 

/ 

2 

26H  SHEAR  WAVE  COR. 

K2 

,  F8 . 6 , 1 5H 

IN**2 

/ 

3 

26H  WAVE  OF  SHEAR  DISC. 

Z 

,F10. 2, 17H 

IN/SEC  (Z-K2*C2)/) 

904 

F0RMAT(24H0R-0  MOMENT  PULSE 

MRO- 

,F8.0,16H 

IN-LB 

/ 

1 

24H  R-0  SHEAR  PULSE 

QRO- 

,  F8. 0, 16H 

LB 

/ 

4 

24H  MRO  PULSE  DURATION 

Tl- 

,  E 1 2 . 5 , 1 4H 

SEC 

/ 

6 

24H  QRO  PULSE  DURATION 

T3- 

,E12. 5 , 14H 

SEC 

/) 

905 

F ORMAT ( 2 6H0TTME  INCREMENT 

DT 

,  E8. 2, 16H 

SEC 

/ 

1 

26H  CELL  LENGTH 

DR 

,F8.5,16H 

IN 

/ 

2 

26H  TIME  SCALE  FACTOR 

TS 

,E8. 2, 16H 

SEC 

/ 

3 

26H  SHEAR  SCALE  FACTOR 

QRS 

,E8. 2 , 16H 

LB 

/ 

4 

26H  VEL  SCALE  FACTOR  VELS 

,E8. 2, 16H 

IN/SEC 

/ 

5 

26H  ANG  VEL  SCALE  FAC  OMGS 

,E8. 2, 16H 

RAD/ SEC 

/ 

6 

26H  MOM.  SCALF.  FACTOR 

MRS 

,E8.2,16H 

LB 

/ 

7 

26H  MOM.  SCALE  FACTOR  MTHS 

, E8. 2, 16H 

LB 

/ 

8 

26H  DEF  SCALE  FACTOR 

WS 

, E8. 2, 16H 

IN 

/) 

906  F ORMAT ( 2  6H0N0 .  OF  CELLS 

IM 

, 14, 15H 

/ 

1 

26H  TOT.  TIME  INTERVALS 

KM 

,  14 , 1 5H 

/ 

2 

26H  NO.  COLS.  PRINTED 

IPR 

,14, 

////) 

907  FORMAT(I4,F8.4) 

908  F ORM AT (  4  6H 1  TIME  PULSE 

1909  F ORMAT( 4 6H 1 TIME  PULSE 

1910  F0RMAT(46H1TIME  PULSE 
911  F0RMAT(46H1TIME  PULSE 
1912  F ORMAT ( 4 6H l TIME  PULSE 
913  FORMAT(lHl) 

4913  F0RMAT(46H1TIME  PULSE 
500  REWIND  21 

DO  127  K-l.KM 
READ  (21)  TT,MMR0,QQ0,(WW(I),I-1,IM0) 
XB(K)»TT 
YB(K)-WW(1) 

YB1(K)-WW<11) 

YB2(K)-WW(21) 

127  CONTINUE 

CALL  I NITQ( A , B , C , D , 600 ) 

CALL  STSWQ(4662,30) 

READ (5 ,800)  (LABEL(J, 1) , J-l ,40) 
READ(5,800)  (LABEL(J,2) ,J-l ,40) 

SET  PLOTIT  ARGUMENT  VALUES. 

N-KM 

LOCX-O 

LOGY-O 


RADIAL  MOMENT  PROPAGATION  /) 

SHEAR  PROPAGATION  /) 

LINEAR  VELOCITY  PROPAGATION  /) 

DEFLECTION  PROPAGATION  /) 

ANGULAR  VELOCITY  PROPAGATION  /) 

TANGENTIAL  MOMENT  PROPAGATION  /) 


i  *  . 
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XAXIS-8.0 

YAXIS-5.0 

XMIN-0.0 

NXDEC-5 

XINC-20.0 

YMIN— 10.0 

NYDEC-5 

YINC-4.0 

XSMIN-0.0 

XSMAX-100.0 

YSMIN— 10.0 

YSMAX-10.0 

NDECX-1 

NDECY-1 

HT-0. 15 

HTS-0.15 

ICON-O 

CALL  PLOTIT(YB , N , LOGX , LOGY , XAXI S , YAXI S , XMIN , XI NC , NXDEC , 

1  YMIN, YINC,  NYDEC, LABEL, XSMIN,XSMAX, YSMIN, YSMAX.XB, SYMA, ICON, HT, 
1  NDECX , NDECY , HTS ) 

DO  130  I-l,N 
130  YB(I)-YB1(I) 

ICON-O 

CALL  SAMEP(YB , N, XB , SYMB , ICON) 

DO  71  I-l.N 
71  YB(I)-YB2(I) 

ICON-O 

CALL  SAMEP(YB , N,XB ,SYMC, ICON) 

CALL  EPLOT 
REWIND  12 
DO  128  K-l.KM 

READ  (12)  TT,MMRO,QQO,(QQ(I ) , I-l , IMO) 

XB(K)-TT 

YB(K)-QQ(1) 

YBKK)-QQ(ll) 

YB2(K)-QQ(21) 

128  CONTINUE 

CALL  INITQ(A,B,C,D,600) 

CALL  STSWQ(4662 , 31 ) 

READ(5,800)  ( LABEL( J , 1 ) , J-l , 40) 

READ(5,800)  (LABEL(J,2) , J-l, 40) 

C  SET  PLOTIT  ARGUMENT  VALUES. 

N-KM 

LOGX-O 

LOGY-O 

XAXIS-8.0 

YAXIS-5 . 0 

XMIN-0.0 

NXDEC-5 

XINC-20.0 

YMIN— 1 0.0 

NYDEC-5 

YINC-4.0 

XSMIN-0.0 
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XSMAX-100.0 

YSMIN— 10.0 

YSMAX-10.0 

NDECX-1 

NDECY-1 

HT-0.15 

HTS-0.15 

ICON-O 

CALL  PLOTIT(YB , N, LOGX , LOGY , XAXIS , YAXIS , XMIN, XI NC , NXDEC , 

1  YMIN,YINC,NYDEC, LABEL, XSMIN.XSMAX, YSMIN, YSMAX.XB.SYMA, ICON, HT, 
1  NDECX.NDECY.HTS) 

DO  72  I-1,N 

72  YB(I)-YB1(I) 

ICON-O 

CALL  SAMEP(YB , N, XB , SYMB , ICON) 

DO  73  I-l.N 

73  YB(I)-YB2(I) 

ICON-0 

CALL  SAMEP(YB , N, XB, SYMC , ICON) 

CALL  EPLOT 
REWIND  11 
DO  129  K-l.KM 

READ  (11)  TT,MMRO,QQO,(MMR(I) ,1-1 ,IM0) 

XB(K)-TT 
YB(K)-MMR( 1 ) 

YBl(K)-MMRdl) 

YB2 (K )«MMR(2 1 ) 

129  CONTINUE 

CALL  INITQ(A,B, C,D, 600) 

CALL  STSWQ(4662 , 32) 

READ(5 , 800)  (LABEL(J , 1 ) , J-l , 40) 

READ(5 , 800)  (LABEL(J , 2) , J-l , 40) 

C  SET  PLOTIT  ARGUMENT  VALUES. 

N-KM 

LOGX-O 

LOGY-O 

XAXIS-8.0 

YAXIS-5 . 0 

XMIN-0 . 0 

NXDEC -5 

XINC-20.0 

YMIN— 10.0 

NYDEC-5 

YINC-4.0 

XSMIN-0.0 

XSMAX-100.0 

YSMIN— 10.0 

YSMAX-10.0 

NDECX-l 

NDECY-l 

HT-0.15 

HTS-0.15 

ICON-O 

CALL  PLOTIT(YB , N, LOGX, LOGY , XAXIS , YAXIS , XMIN, XI NC , NXDEC, 
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1  YMIN.YINC,  NYDEC ,  LABEL,  XSMIN ,  XSMAX .YSMIN ,  YSMAX ,  XB ,  SYMA ,  ICON,  HT , 
l  NDECX, NDECY.HTS) 

DO  74  I-l,N 
74  YB(I)-YBt(I) 

ICON-O 

CALL  SAMEP(YB , N, XB , SYMB , ICON) 

DO  75  I-l.N 
75  YB(I)-YB2(I) 

ICON-O 

CALL  SAMEP (YB , N, XB , SYMC , ICON) 

CALL  EPLOT 
REWIND  21 
DO  131  K-l.KM 

READ  (21)  TT.MMRO, QQO , ( WW(I) ,1-1 ,  IHO) 

IF(K. NE.KKM. AND.K. NE .KKM1 . AND.K. NE .KKM2)  GO  TO  131 

DO  231  I-l.IHO 

XB(I)-R(I) 

IF(K.EQ.KKM)  YB(I)-WW(I) 

IF(K.EQ.KKMl)  YB1(I)-WW(I) 

IF(K.EQ.KKM2)  YB2(I)-WW(I) 

231  CONTINUE 
131  CONTINUE 

CALL  INITQ(A,B,C,D,600) 

CALL  STSWQ(4662 , 33) 

READ(5, 800)  (LABEL(J , 1) ,J-1 ,40) 

READ(5 ,800)  (LABEL(J ,2) ,J-1 ,40) 

C  SET  PLOTIT  ARGUMENT  VALUES. 

N-IMO 

LOGX-O 

LOGY-O 

XAXIS-8.0 

YAXIS-5.0 

XMIN-0.0 

NXDEC-5 

XINC-O. 2*DS 

YMIN— 10.0 

NYDEC-5 

YINC-4.0 

XSMIN-0 . 0 

XSMAX-DS 

YSMIN— 10.0 

YSMAX-10.0 

NDECX-1 

NDECY-l 

HT-0.15 

HTS-0.15 

ICON-O 

CALL  PLOTIT(YB , N, LOGX , LOGY , XAXI S , YAXI S , XMIN , XI NC , NXDEC , 

1  YMIN.YINC , NYDEC , LABEL, XSMIN, XSMAX, YSMIN, YSMAX, XB, SYMA, ICON, HT, 
1  NDECX,. NDECY.HTS) 

DO  76  I-l.N 
76  YB(I)-YB1(I) 

ICON-O 

CALL  SAMEP(YB,N,XB, SYMB, ICON) 
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DO  77  I«l,N 
77  YB(I)-YB2(I) 

ICON-O 

CALL  SAMEP(YB,  N.XB.SYMC, ICON) 

CALL  EPLOT 
REWIND  12 
DO  132  K-l,KM 

READ  (12)  TT.MMRO, QQO ,  (QQ(I ) , 1-1 , IMO) 

IF(K. NE.KKM.AND.K. NE.KKM1.AND.K. NE.KKM2)  GO  TO  132 

DO  232  1*1 , IMO 

XB(I)-R(I) 

IF(K.EQ.KKM)  YB(I)-QQ(I) 

IF(K.EQ.KKMl)  YB1(I)-QQ(I) 

IF(K.EQ.KKM2)  YB2(I)*QQ(I) 

232  CONTINUE 
132  CONTINUE 

CALL  INITQ(A,B,C, D.600) 

CALL  STSWQ(4662 , 34) 

READ ( 5 , 800 )  (LABEL(J,1),J-1,40) 

READ(5 , 800)  (LABEL(J , 2) , J-l ,40) 

C  SET  PLOTIT  ARGUMENT  VALUES. 

N-IMO 

LOGX-O 

LOGY-O 

XAXIS-8.0 

YAXIS-5.0 

XMIN-0.0 

NXDEC-5 

XINC-O. 2*DS 

YMIN— 10.0 

NYDEC-5 

YINC-4.0 

XSMIN-0.0 

XSMAX-DS 

YSMIN—10.0 

YSMAX-10.0 

NDECX-1 

NDECY-1 

HT-0.15 

HTS-0.15 

ICON-O 

CALL  PLOTIT(YB , N, LOGX , LOGY , XAXI S , YAXI S , XMIN , XI NC , NXDEC , 

1  YMIN, YINC.NYDEC, LABEL, XSMIN,XSMAX,YSMIN,YSMAX,XB,SYMA, ICON, HT, 
1  NDECX , NDECY , HTS ) 

DO  78  1*1 , N 
78  YB(I)-YB1(I) 

ICON-O 

CALL  SAMEP ( YB ,  N,  XB , SYMB , ICON) 

DO  79  1*1 , N 
79  YB(I)-YB2(I) 

ICON-O 

CALL  SAMEP(YB , N, XB , SYMC , ICON) 

CALL  EPLOT 
REWIND  11 
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DO  133  K-l.KM 

READ  (11)  TT,MMRO,QQO,(MMR(I) ,1*1 ,IMO) 

IF(K. NE.  KKM.  AND. K.  NE . KKM1 . AND. X. NE . KKM2 )  GO  TO  133 

DO  233  I-l.IMO 

XB(I)-R(I) 

IF(K.EQ.KKM)  YB(I)-MMR(I) 

IF(K.EQ.KKMl)  YB1(I )-MMR(I ) 

IF(K.EQ.KKM2)  YB2(I )-MMR(I ) 

233  CONTINUE 
133  CONTINUE 

CALL  INITQ(A, B,C,D, 600) 

CALL  STSWQ(4662,35) 

READ(5,800)  (LABEL(J,l),J-i,40) 

READ(5 , 800)  (LABEL(J , 2) , J-l , 40) 

C  SET  PLOTIT  ARGUMENT  VALUES. 

N-IMO 

LOGX-0 

LOCY-O 

XAXIS-8.0 

YAXIS-5.0 

XMIN-0.0 

NXDEC-5 

XINC-O. 2*DS 

YMIN— 10.0 

NYDEC-5 

YINC-4.0 

XSMIN -0.0 

XSMAX-DS 

YSMIN— 10.0 

YSMAX-10.0 

NDECX-1 

NDECY-l 

HT-O. 15 

HTS-0.15 

ICON-O 

CALL  PL0TIT(YB , N, LOGX , LOGY , XAXIS , YAXIS , XMIN, XI NC , NXDEC , 

1  YMIN, YINC.NYDEC,  LABEL, XSMIN, XSMAX,  YSMIN, YSMAX.XB.SYMA,  ICON, HT, 
1  NDECX , NDECY , HTS ) 

DO  81  I-l.N 
81  YB(I)-YBKI) 

ICON-0 

CALL  SAMEP  (YB ,  N,  XB ,  SYMB ,  ICON) 

DO  82  I-l.N 
82  YB(I)-YB2(I) 

ICON-O 

CALL  SAMEP  (YB ,  N,  XB ,  SYMC ,  ICON) 

CALL  EPLOT 
YB(l)-0.0 
YBl(l)-0.0 
YB2(1 )-0.0 
REWIND  21 
DO  136  K-l.KM 

READ  (21)  TT,MMR0,QQ0,(WW(I),I-1 , IMO) 

IF(K. NE.XXM.AND.K. NE.KKM1. AND.K. NE.KKM2)  GO  TO  136 
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DO  236  1-1 ,IM 
XB(I)-R(I) 

IF(K. EQ. KKM)  YB( I+l )-(WW( I+l )-WW( I ) ) /DR 
IF(K.SQ.KKMl)  YB1 (I+l )-(WW( I+l )-WW(I ) )/DR 
IF(K.EQ.KKM2)  YB2(I+1 )-(WW(I+l)-WW(I))/DR 
236  CONTINUE 
136  CONTINUE 

CALL  INITQ(A,B,C,D,600) 

CALL  STSWQ(4662 , 36) 

READ(5 , 800)  (LABEL(J , 1 ) , J-l, 40) 

READ(5 , 800)  (LABEL(J,2) , J-l ,40) 

C  SET  PLOTIT  ARGUMENT  VALUES. 

N-IMO 

LOGX-O 

LOGY-O 

XAXIS-8.0 

YAXIS-5 . 0 

XMIN-0.0 

NXDEC-5 

XINC-O. 2*DS 

YMIN— 10.0 

NYDEC-5 

YINC-4.0 

XSMIN-0.0 

XSMAX-DS 

YSMIN— 10.0 

YSMAX-10.0 

NDECX-1 

NDECY-1 

HT-O. 15 

HTS-0.15 

ICON-O 

CALL  PLOTIT(YB , N, LOGX , LOGY , XAXIS , YAXIS , XMIN, XI NC , NXDEC , 

1  YMIN,  YINC ,  NYDEC ,  LABEL,  XSMIN,  XSMAX,  YSMIN,  YSMAX,  XB ,  SYMA ,  ICON,  HT , 

1  NDECX, NDECY , HTS ) 

DO  84  I-l.N 
84  YB(I)-YBl(I) 

ICON-O 

CALL  SAMEP(YB ,  N, XB , SYMB , ICON) 

DO  85  I-l , N 
85  YB(I)-YB2(I) 

ICON-O 

CALL  SAMEP(YB , N, XB , SYMC , ICON) 

CALL  EPLOT 

STOP 

END 

SUBROUTINE  PLATE 

DIMENSION  R(250)  ,QR(250)  ,VEL(250)  ,W(250)  ,0MEGA(250)  ,QQ(250) , 
1VEU7E(250),0M0M(250),WW(250),AIMR(250),AIMTH(250),U(250),AIQR(250) 
1,QQS(250) 

REAL  MR(250),MTH(250),NU,MRO,MRS,MTHS,MMR(250),MMTH(250), 

1MMR0 , K2 , LIQR(2  50 ) 

LOGICAL* 1  IDENT(80) ,FMT(80) ,NAME(10) 

COMMON  T,R,G,EM,IMO,H,DI,RO,RL,MR,MTH,QR,CP,C2,K2,VEL,M,RHO, 
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1 DE  N, NU , OMEGA , MRO , QRO , T1 , T3 , DS  TR  N , DALPH , DT , DR , TS , MRS , 

2MTHS , QRS , V E LS , OMGS , WS , TT , QQ , MMR , MMTH , V E LV R , OMCM , IM , KM , I PR , 
3WV , MMRO , QQO , QQS , 2 , U , TAU1 , TAU2 , ESR , NAME , FMT , IDE  NT 
PRINT  921 
921  FORMAT(lHl) 

C  INITAL  STATE  OF  PLATE 
DO  119  I-l.IMO 
AIMR(I)-0.0 
AIMTH(I)-0.0 
AIQR(I)-0.0 
LIQR(I)-0.0 
MR(I )-0 . 0 
MTH(I)»0.0 
QR(I)-0.0 
VEL(I )-0 . 0 
W(I)-0.0 
0MEGA(I)-0.0 
U(I)»0.0 
119  CONTINUE 
DTD-DR/CP 
DTR-DR/Z 
T-0.0 
TD-0.0 
TL-0 . 0 
TLB-0.0 
TR-0.0 

DO  101  K-l.KM 
TT-T/ TS 

C  ***INPUT  AT  R“RO,  MOMENT  BOUNNDARY  CONDITION*** 

IF(T-Tl)  110,111,111  HO  MR(IMO)-MRO 
GO  TO  1492  111  MR(IM0)-0.0 

1492  CONTINUE 

C  ***INFUT  AT  R-RO,  SHEAR  BOUNDARY  CONDITION*** 

IF(T-T3)  113,114,114 

113  QR(1)-QR0/(2*3.1415*R(1)) 

GO  TO  115 

114  QR(l)-0,00 

115  CONTINUE 

C  ***FOUNDATION  REACTION*** 

DO  116  I-l.IMO 
QQS  (I)— ESR*W(I) 

116  CONTINUE 

C  SCALING  FOR  PRINTOUT 

IZZ-IMO 

DO  103  I-l.IZZ 
MMR(I)-MR(I)/MRS 
MMTH(I)-MTH(I)/MTHS 
QQ(I)»QR(I)/QRS 
VELVE(I)-VEL(I)/VELS 
OMOM( I )-OMEGA( I ) /OMGS 
WW(I)-W(I)/WS 
103  CONTINUE 

WRITE ( 1 1 )TT , MMRO , QQO , ( MMR( I ) , I -l , IMO ) 

WRITE ( l 2 )TT , MMRO , QQO , ( QQ( I ) , 1-1 , IMO) 
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WRITE(13)TT,MMR0,QQ0,(MMTH(I),I-1 , IMO) 

WRITE(20)TT,MMR0, QQO ,  (7  ELVE  (I  )  ,  1-1 ,IMO) 

WRITE/ 21 )TT, MMRO , OOO , ( WW( I ) ,I«l , IMO) 

WRITE(22)TT,MMRO,QQO,(OMOM(I) ,1-1 , IMO) 

12  IF(ABS(TD-(T+0T))-0. IE-10)  19,19,89 
89  IF(TD-TR)  15,15,16 

16  IF(TD-M)TD-TR-DTR)  15,15,13 
C  THE  CP  WAVE 

15  DO  108  1-2, IMO 

AIMR(I)-AIMR(I)+(MR(I)*(R(I-1 )+DR)-MR(I-l )*R(I-1 ))*DTD 
AIMTH(I)-AIMTH(I)+MTH(I-1 )*DR*DTD 

AIQR(I)-AIQR(I)+(QR(I)*(R(I-1 )+DR)*DR+QR(I-l )*R(I-1 )*DR)*0.5 
l *(TD+DTD-TL) 

LIQR/I-1 )-LIQR(I-l )+(QR(I)*(R(I-l )+DR)-QR(I-l )*R(I-1 ) 

1  )*(TD+DTD-TL) 

IF(TAU1-1 .0E+40)  930,928,928 
930  DOMG-(-l . 0)*OMEGA( I )*(DTD/TAU1 ) 

OMEGA(I )-OMFCA( I )+(D0MG/2 .0) 

C  IMPULSE  AND  MOMENTUM  ACROSS  CP. 

928  B-((RHO*H)/24.0)*((H**2))*((2.0*R(I)*DR)+(DR**2)) 

DOM-(AIMR(I)-AIMTH(I)-AIQR(I))/B 
OMEGA( I )-OMEGA(I )+DOM 
OMEGA ( 1 )«0 . 0 

IF( TAU1-1 . 0E+40)  9113,815,315 
C  DAMPING,  SECOND  HALF 

9113  D0MG-/-1 . 0)*0MEGA( I )*(DTD/TAUl ) 

0MEGA( I )-0MEGA( I )+(D0MC/2 . 0) 

C  INITIALIZE  ANGULAR  IMPULSE  AFTER  USE 

815  AIMR(I )-0. 0 
AIMTH(  I  )-0 . 0 
AIQR(I)-0. 0 
108  CONTINUE 
TD-TD+DTD 
TLB-TL 
TL-TD 

DO  107  I-l.IM 
C  ANGULAR  STRAIN 

DALPH-( (OMEGA( 1+1 )-OMEGA( I ))*DTD ) /DR 
C  CONSTITUTIVE  EQUATIONS  FOR  MOMMENT 

IF(I.EQ.l)  DMR-D I * ( D A LPH ) 

IF(I.NE.l)  DMR-DI*(DALPH+((NU)*U(I-1 ))) 

DM  TH-D I  *  (  U  ( I H  (  (  NU  )  *D  A  LPH  )  ) 

MR(I)-MR(I)+DMR 
IF(I.EQ. 1)  MTH(I)-MTH(I)+DMTH 
IF(I.NE.l)  MTH(I-1 )“MTH(I-1 )+DMTH 
U ( I )- ( OMEGA( I ) *DTD ) / R( I ) 
QR(I+l)-QR(I+l)+(K2**2)*G*H*OMEGA(I+l)*DTD 
107  CONTINUE 

GO  TO  12 

13  DO  104  I-l.IM 

LIQR(I)-LIQR(I)+(QR(I+1)*(R(I)+DR)-QR(I)*R(I))*(TR+DTR-TL)+ 
l (0. 5*(QQS(I+l )+QQS(I ) )*(R(I )+DR/2 . 0)*DR)*DTR 
AIQR(I+1 )-AIQR(I+l)+(QR(I+l)*(R(I)+DR)*DR+QR(I)*R(I)*DR)*0.5*(TR+ 
1DTR-TL) 
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IF(TAU2-1 .0E+40)  925,931,931 
C  DAMPING,  FIRST  HALF 
925  DVELG«(-1 .0)*/EL(T)*(DTR/TAU2) 

VEL(I)-VEL(I)+(DVELG/2.0) 

C  IMPULSE  AND  MOMENTUM  ACROSS  K2*C2 
931  A-<R(I)*DR)+((0.5)*DR*DR) 

WEE-LIQR(I)/(RHO*H*A) 

VEL<T)^/EL(I)+DVEE 

VEL(IM+1)«0.0 

IF(TAU2-1 . 0E+40)  927,818,818 
C  DAMPING,  SECOND  HALF 
927  EWELG-(-1.0)*VEL(I)*(DTR/TAU2) 

VEL(I)*/EL(I)+(DVF.LG/2.0) 

818  LI OR (I )-Q.  0 
104  CONTINUE 
TR-TR+DTR 
TLB-TL 
TL-TR 

DO  102  I-l.IM 

DSTRN-((VEL(I+1  )-VEL(I)  )*DTR)/DR 
DOR»(K2**2 )*G*H*(DSTRN) 

QR(I+1)-QR(I+1 )+DQR 
W(I)-W(I)+0/EL(I)*DTR) 

102  CONTINUE 
GO  TO  12 
19  CONTINUE 
T-T+DT 

101  CONTINUE 
ENDFILE  11 
ENDFILE  12 
ENDFILE  13 
ENDFILE  20 
ENDFILE  21 
ENDFILE  22 
500  RETURN 
END 

//DATA. FTilFOOl  DD  UNIT-SYSDA, SPACE-(TRK, ( 5, 5) ,RLSE) 

//DATA.FT12F001  DD  UNIT-SYSDA, SPACE-(TRK, ( 5 , 5) ,RLSE) 

//DATA. FT13F001  DD  UNIT-SYSDA, SPACE-(TRK, (5, 5) ,RLSE) 

/ /DATA. FT20F 001  DD  UNIT-SYSDA, SPACE-(TRK, (5, 5) ,RLSE) 

//DATA.FT21F001  DD  UNIT-SYSDA, SPACE-(TRK, (5 , 5) ,RLSE) 

//DATA.FT22FOOI  DD  UNIT-SYSDA, SPACE-(TRK, (5 , 5) ,RLSE) 

//DATA. INPUT  DD  * 

CONDITION  1  FOR  PLATE  ON  ELASTIC  FOUNDATION 
(IHS,F6.2,1X,F5.1,1X,F5.1,1X,22F5.2) 

CONCRETE  0.0868  30.00E+5  0.20  15.00  0.20  0.20 

4.00 

0 . OOOE+O  -100.0E+0  4 . OOOOE+O  4.0000E+0  3.0000E-5  3.0000E-5 

l.OOOOE-5  10. OOOE+O  4.0000E-1  5.0000E-6  5.00E+0  l.OOOOE-K)  5.0000E-K) 

100  22  100  100  50 

1.20E+5  15.00E+0  1 

TIME  IN  1.0*E-5  SEC 
PLATE  DEFLECTION  IN  5*E-6  INCHES 
TIME  IN  l.0*E-5  SEC 
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PLAT1!:  SHEAR  IN  10.*LB. 

TIME  IN  1 . 0*E-5  SEC 
PLATE  MOMENT  IN  5.*LB-IN. 

DISTANCE  X  IN  INCHES 

PLATE  DEFLECTION  IN  5*E-6  INCHES 

DISTANCE  X  IN  INCHES 

PLATE  SHEAR  IN  10.*LB 

DISTANCE  X  IN  INCHES 

PLATE  MOMENT  IN  5.*LB-IN. 

DISTANCE  X  IN  INCHES 
PLATE  SLOPE  IN  5*E-6 
//DATA. ET30F001  DD  DSN-MEN. P11850. A3B 
//  VOL-REF-MEN.P11850.A3B 

//  SPACE-(TRK, ( 1 ,2) ,RLSE) 

//  DCB-(RECEM-FB , LRECL-30 

//DATA.FT31F001  DD  DSN-MEN. PI 1850. A3B 
//  VOL-REF-MEN. P11850. A3B 

//  SPACE-(TRK, (1,2), RLSE) 

//  DCB-(RECFM-FB, LRECL-80 

/ /DATA. FT32F001  DD  DSN-MEN. PI 1850. A3B 
//  VOL-REF-MEN. P11850.A3B 

//  SPACE-(TRK,( 1,2), RLSE) 

//  DCB-(RECFM-FB, LRECL-80 

/ /DATA. FT33F001  DD  DSN-MEN. P11850. A3B 
//  V0L-REF-MEN.P11850.A3B 

//  SPACE-(TRK,( 1,2), RLSE) 

//  DCB-(RECFM-FB, LRECL-80 

/ /DATA . FT34F001  DD  DSN-MEN. PI 1850 . A3B 
//  VOL-REF-MEN. PI 1 850. A3B 

//  SPACE-(TRK,( 1,2), RLSE) 

//  DCB-(RECFM-FB, LRECL-80 

/ /DATA. FT35F001  DD  DSN-MEN. PI 1 850. A3B 
//  V  OL-REF -MEN. PI 1 8  50 . A3B 

//  SPACE- ( TRK, ( 1,2) ,RLSE) 

//  DCB-(RECFM-FB, LRECL-80 

/ /DATA. FT36F001  DD  DSN-MEN. P11850.A3B 
//  VOL-REF-MEN. P11350.A3B 

//  SPACE-(TRK, (1,2), RLSE) 

//  DCB-(RECFM-FB, LRECL-80 


.PLOT, 

.  LIB, DISP-( NEW, KEEP), 

,  BLKSIZE-1 2960) 

.  PL0T1 , 

.  LIB, DISP-( NEW, KEEP) , 

,  BLKSIZE-12980) 

.  PL0T2 , 

.  LIB, DISP-( NEW, KEEP), 

,  BLKSIZE-1 2960) 

.PL0T3 , 

.  LIB, DISP-(NEW, KEEP), 

,*  BLKSIZE-1 2960) 

.  PLOTA , 

. LIB, DISP-( NEW, KEEP), 

,  BLKSIZE-1 2960) 

. PLOT 5 , 

. LIB, DISP-(NEW, KEEP) , 

, BLKSIZE-1 2960) 

. PL0T6 , 

. LIB, DISP-(NEW, KEEP), 
,  BLKSIZE-1 2960) 


